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(57) Abstract 



Apparatus and method are provided for combining, in real time, weather radar data in digital format from multiple radar sites into a 
mosaic covering a regional or national area (e.g., the continental United States). Such apparatus and method is designed to be implementable 
as computer software and to run on a general purpose computer such as a personal computer. Such apparatus and method preferably utilizes 
a database of "lookup tables" which are used to project individual radar data bins directly to a grid of the desired coordinate system for 
that mosaic. Using these lookup tables, each of a plurality of radar data bins is converted to a corresponding grid location or box in a 
mosaic image. While each radar data bin is mapped to a grid box for the mosaic, not every such grid box within the coverage area of the 
radar may receive such a data bin. Any resulting holes in the produced mosaic are then identified using yet another grid whose locations 
correspond to those of the mosaic but which instead provide information regarding the data then appearing in corresponding locations in 
that mosaic. Each such hole is then filed with data from the locationally closest data bin. 
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APPARATUS AND METHOD FOR CQNBTRPCTTttfl 

ft ttflgAIS OF DATA 

A portion of the disclosure of this patent document 
contains material which is subject to copyright protection. 
The copyright owner has no objection to the facsimile 
reproduction by anyone of the patent document or the patent 
disclosure, as it appears in the U. S. Patent & Trademark 
Office patent files or records, but otherwise reserves all 
copyright rights whatsoever. 

FIELD OF THE INVENTION 
The present invention relates to display of radar 
images, and more particularly to producing such a display 
by combining a plurality of radar images into a mosaic of 
those images. The present invention also relates to 
electrical computers and data processing systems, and more 
particularly to applications of same to earth science such 
as weather. The present invention further relates to 
information processing system organization, and more 
particularly to applications using artificial intelligence 
with details of the artificial intelligence system such as 
for earth sciences such as weather. The present invention 
still further relates to communications, directive radio 
wave systems and devices (e.g. radar or radio navigation) , 
and more particularly to weather radar, plural radar, 
display circuits for image production or scan conversion 
for same, and directive radio wave systems and devices for 
position indicating (e.g. triangulation) for storm 
location. 

BACKGROUND Of THE INVENTION 
There exist applications in which there is a need to 
combine several radar images into a singl mosaic image 
covering a geographical ar a. Th pr sent invention 
fulfills that need. 
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There is a large number of weather radar sit s 
distributed across the U.S.A. There is a national weath r 
radar network consisting of the new NEXRAD WSR-88D weather 
surveillance Doppler radars and the preexisting WSR-57 and 
WSR-74 non-Doppler weather radars. Each such site can 
essentially operate independently, but the ranges of thes 
sites can overlap. It is necessary to assemble the data 
from these sites into a single mosaic showing weather 
conditions in a particular region or across the country. 
The present invention fulfills this need. 

A mosaic of radar data can provide location and 
intensity of precipitation and severe weather throughout a 
large geographic area and is, therefore, useful to 
meteorologists and other persons responsible for monitoring 
and forecasting weather conditions. Creation of a mosaic 
of radar data usually requires conversion of data points 
from a local geographic coordinate system (wherein the 
positions there are defined relative to an individual radar 
site) to the geographic coordinate system of the mosaic. 
Conversion of data points from one coordinate system to 
another coordinate system (reprojection) is a 
computationally intensive task, but is necessary to ensure 
accuracy. Radar data from even one site, however, can 
consist of over 100,000 individual data values, each of 
which must be reprojected. A mosaic could be made up of 
data from dozens of such radar sites. There is therefore 
a need to accurately combine all of that data from all 
desired sites fast enough to provide a current depiction of 
the weather. The present invention fulfills this need. 

SUMMARY OF THE INVENTION 

Accordingly, an object of the present invention is to 
provide apparatus and method for combining data from a 
plurality of sites into a single image covering a desired 
geographical area. 

Another object of the present invention is to pr vid 
apparatus and method capable of combining, in real time, 
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radar data (such as weather radar data) in digital £ rmat 
from multiple radar sites into a mosaic covering a 
geographical area. 

A further object of the present invention is to 
provide apparatus and method for combining data from a 
plurality of sources into a single image that can be run on 
a general purpose computer such as a personal computer. 

Still another object of the present invention is to 
provide apparatus and method capable of producing a mosaic 
of radar data that can provide location and intensity of 
precipitation and other meteorological phenomena throughout 
a large geographic area. 

A still further object of the present invention is to 
provide apparatus and method capable of accurately 
combining a substantial amount of data from a plurality of 
desired sites fast enough to provide a current depiction of 
the weather or other phenomena. 

Yet another object of the present invention is to 
provide apparatus and method capable of accurately yet 
efficiently building mosaics of radar data. 

Briefly, these and other objects of the present 
invention are accomplished by apparatus and method for 
converting each of a plurality of radar data bins to a 
corresponding grid location or box in a mosaic image. 
However, while each radar data bin is mapped to a grid box 
for the mosaic, not every such grid box within the coverage 
area of the radar may receive a bin. This may occur 
because there may be more grid boxes available than there 
are radar bins, or because the radar data bins do not 
overlap perfectly. After the mapping from bins to boxes is 
completed, the inventive apparatus and method then goes 
back and fills those holes in the following manner. Bins 
of radar data are mapped to a predetermined grid for a 
portion of the resulting display. Each time that a grid 
element f that grid receives a bin, then a corresponding 
element in th temporary array is "set". The temporary 
array is then ch eked for any "unset" or blank el ments. 
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If an array el m nt is found without a bin, then the 
locationally closest data bin is found. That closest data 
bin and its corresponding grid location are used to fill 
the hole in the mosaic grid. 

Other objects , advantages and novel features of th 
invention will become apparent from the following detail d 
description of the invention when considered in conjunct! n 
with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWTKGS 

In the drawings , 

Fig. 1 is a diagrammatic representation of a tw 
dimensional array for the display of a single site radar 
product; 

Fig. 2 is a diagrammatic representation of one exampl 
of mapping of individual weather radar product data into a 
larger mosaic grid; 

Fig. 3 is a diagrammatic representation of an overview 
of a process according to the present invention to 
indirectly map individual radar data values to a mosaic 
grid according to the present invention; 

Fig. 4 is a diagrammatic representation of an overview 
of a process to directly map individual radar data valu s 
to a mosaic grid according to a preferred embodiment of the 
present invention; 

Fig. 5 illustrates the format of a lookup table that 
can be utilized in the method of Fig. 4; 

Fig. 6 shows one illustrative example of a display 
resulting from radar data mapping that has holes present 
where the source data grids and the destination data grid 
do not match up; 

Figs. 7A and 7B together provide a diagrammatic 
representation of one embodiment of a method of filling any 
one of the holes of Fig. 6 according to the present 
invention; 
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Figs. 8A, 8B, 8C, 8D and 8E tog ther sh v a f 1 wchart 
for one method of building the lookup table of Figs. 4 and 
5 ; and 

Figs. 9A, 9B, 9C, 9D and 9E together show a flowchart 
for one method of building a mosaic of radar data according 
to the present invention using the lookup tables of Figs. 
4, 5 and 8 A through 8£. 

DETAILED DESCRIPTION 

Referring now to the drawings, wherein like reference 
characters designate like or corresponding parts throughout 
the several views, there is shown in Fig. 1 an array 10 
made up of a plurality of elements or pixels, of which one 
is identified by reference numeral 12. Each array element 
( such as element 12 ) contains a single value of radar 
return. Each such element represents a finite area on a 
horizontal plane essentially tangent to the earth's 
surface. Array 10 contains one example of a decoded single 
site radar data product 14. In array 10, product 
representation 14 is bounded by circle 16 representing the 
extent of weather radar data return within the product. 
Within product 14, an individual radar data bin location 
(such as location 18) is identified by its x, y coordinates 
in the array, with the bin latitude and longitude being 
computed based on its distance relative to the radar site. 
Central bin 20 represents the radar site location. For a 
fixed radar site, the latitude and longitude of the radar 
site is known, so that referencing an individual bin with 
respect to location of the central bin will provide precise 
position information for that individual bin also. The 
single site radar data is received in digital format which 
allows the data values to be assigned to a computer 
program's memory in a two-dimensional array 10. Each array 
element (such as element 12) is a data "bin" which 
represents an atm spheric parameter (e.g. pr cipitation 
intensity) within a relatively small finit area such as 1 
kil meter by 1 kilometer. Thus, the coordinates f any one 
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bin (such as bin 18) can provide the g ographic x,y 
distance of that bin from the radar site. 

However, often more than one radar product must be 
used in a single grid array or visual representation of 
weather conditions in a region. This may be because the 
range for two radar sites may overlap, or because it is 
desired to show, in one display, radar data for a 
geographic area that is sufficiently large to includ a 
plurality of radar sites. In such situations, a mosaic is 
built by taking the individual data bins from the radar 
data for a single site and assigning the values for thos 
bins to the corresponding locations in a large scale (e.g. 
regional or national) grid which preferably also resides in 
a computer memory as a two-dimensional array 22. A simpl 
example of formation of such a mosaic is shown in Fig. 2. 
Fig. 2 shows one example of mapping of individual weath r 
radar product data from arrays 10, 24 and 26, such as from 
different radar sites, into a larger scale mosaic grid 22. 
Each array 10, 24 and 26 here represents radar product from 
a single site A, B, or C respectively. Each such array 10, 
24 and 26 is similar to array 10 of Fig. 1. As shown in 
Fig. 2, the data from array 10 is combined with data from 
other such arrays 24 and 26 to construct a mosaic grid 22 
which covers the geographic area that includes all thre 
such radar sites. Thus, in the example of Fig. 2, array 10 
is mapped onto portion or subarray 28 of grid 22; array 24 
is mapped onto portion or subarray 30 of grid 22; and array 
26 is mapped onto portion or subarray 32 of grid 22. In 
the example of Fig. 2, the coordinate system of the radar 
data from a single radar radar does not necessarily match 
the coordinate system of the mosaic grid 22. That is, th 
bins of arrays 10, 24 and 26 and the grid boxes of grid 22 
are not necessarily the same size, and the respectiv 
coordinate systems are not necessarily oriented the same 
way. For this reas n, the co rdinates of each bin fr m an 
array 10, 24 r 26 must be converted (i.e. r proj cted) t 
the coordinat s of grid 22. 
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On m th d of reprojecting the coordinates of an 
individual data point such as for transfer from array 10 to 
grid 22 according to the present invention is shown in Fig. 
3. Fig. 3 shows an overview of a process according to the 
present invention to indirectly map individual radar data 
values such as from array 10 to a mosaic grid such as grid 
22. In single site radar product array 10, data values are 
assigned to bins. In the method of Fig. 3, each bin is 
converted to latitude and longitude and then that latitude 
and longitude is converted to the corresponding grid 
location. At step 34, each set of bin coordinates x,y is 
converted to distance from radar site 20 on the x and y 
axes. At step 36, this distance in the form of x and y 
values is then converted to geographic latitude and 
longitude. At step 38, these values of latitude and 
longitude are then converted to x, y distance from the grid 
22 center, now using the x and y axes of grid 22. At step 
40, this x, y distance is converted to the corresponding x 
and y coordinate values for grid 22. In this manner, data 
values from an array 10 are mapped to mosaic grid 22. A 
similar procedure can be utilized for mapping from other 
arrays 24, 26 to grid 22. Although this approach is 
accurate, it is not preferred because it is computationally 
intensive (see steps 34, 36, 38 and 40) and therefore would 
not be desirable for real-time operation on a general 
purpose computer such as a personal computer. 

An alternative, more preferred method for mapping from 
a single site grid or product array 10 to a mosaic grid 22 
is illustrated in Fig. 4. The method of Fig. 4 is more 
efficient in real time than that of Fig. 3 because the 
method of Fig. 4 maps the radar data bins directly to the 
grid 22 using prebuilt look-up tables. These look-up 
tables are built using computationally intensive 
algorithms, illustrated in Fig. 8 and shown in the 
c rresponding code for that figure, to compute latitude and 
longitude f each radar data bin and to convert th se 
latitude and longitud valu s to b x locations f grid 22. 
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specific to th coordinate system and spatial r solution of 
the individual product (e.g. product 14) and of the mosaic 
grid (e.g. grid 22). A given bin x, y coordinate can be 
paired with more than one grid x, y coordinate. Likewise, 
a given grid x, y coordinate can be paired with more than 
one bin x, y coordinate. However , there is at least one 
table entry for each bin x, y coordinate, and there is at 
least one entry for each grid x, y coordinate. However, 
the total number of entries likely would not exceed the 
larger of the number of bins or the number of boxes to be 
mapped • 

As discussed above, the array of product data bins 
(e.g. array 10 , 24 or 26) will not necessarily be perfectly 
aligned with mosaic grid 22. As a result, some grid boxes 
will receive more than one bin. For similar reasons, some 
grid boxes, within the radar data range, will not receive 
any bins. In the latter case, "holes" would remain in the 
completed mosaic, such as can be seen in Fig. 6. Fig. 6 
shows one example of radar data mapping having one or more 
holes 52 where the source and destination grids, i.e. the 
radar bins and mosaic grid, do not match up. The following 
additional processing is utilized in the present invention 
to fill such holes 52 in grid 22. This process is 
accomplished in the following manner. In the process to 
build look-up tables as shown in Fig. 8, the location of 
each such "hole" in the grid is identified and the closest 
bin is assigned to the grid box having that hole. This is 
done by using the temporary array of Fig. 7 to find 
unfilled grid boxes within the radar coverage area. Each 
time a bin x, y coordinate is paired with a grid x, y 
position, an element in temporary array 54 is set. Then, 
when processing of all bin x, y locations for a site is 
complete, a "backward calculation" is performed for all 
unset or unfilled grid boxes in temporary array 54 within 
the range of the radar for which a look-up table is then 
being constructed. In such a backward calculation, the 
grid b x x, y locati n is converted to latitud and 
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longitude, and then that latitude and longitude is 
converted to the corresponding bin x, y coordinate to find 
the bin closest in location to that bin x, y coordinate for 
the hole. Each so "back calculated" bin x, y and grid x, 
y pair is added as an entry to the look-up table for that 
radar site. This process ensures that every data bin is 
assigned to the mosaic grid and that every grid box within 
the range of the radar then being considered is assigned a 
data value. 

Fig. 7 shows use of a temporary array 54 to fill a 
hole in grid 22. Each element of temporary array 54 can 
contain any one of three values: a first, initial "blank" 
value (such as -1); a second "unset" value (such as 0) for 
each corresponding grid element or box falling within the 
range of a radar; and a third "set" value (such as +1) 
indicating that a radar bin has been mapped to the 
corresponding grid element. In a subarray (such as 
subarray 32) of grid 22, bins of weather radar data from 
the corresponding product array (e.g. array 26) are mapped 
to grid 22 with the orientation of grid 22 and not that of 
the source array then being mapped. In this mapping 
process, each time that a grid 22 element receives a bin, 
then a corresponding element in temporary array 54 is "set" 
to so indicate. Temporary array 54 is then checked for any 
"unset" or unfilled elements. In temporary array 54, 
circle 56 represents the range of weather data in the 
weather radar product then being mapped. All array 54 
elements within circle 56 are to receive a data bin or 
value from the corresponding product array being mapped to 
the grid 22. If an array element of array 54 is found that 
has not been so "set", then a back calculation is performed 
to find the filled data bin that is closest in location to 
that array element. That data bin x, y and its 
corresponding grid location x, y are added to the look-up 
table to fill the hole that would otherwise result in 
m sale grid 22. 
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To perform the proc ss of Fig. 4, a two-part process 
to build mosaics is needed. The first of these two parts, 
shown in Figs. 8A-8E, can be executed off-line and is 
responsible for the generation of look-up tables which are 
built and stored in a local database 46. The second part, 
which is designed to run in real time and uses the pre- 
built look-up tables to actually build the mosaic (s) , is 
shown in Figs. 9A-9E. 

In Figs. 8A-8E, the process to build site-specific 
look-up tables can be summarized as follows. The following 
process is conducted for each radar site, with a separate 
table being constructed for each radar site. Each element 
in temporary array* 54 is first initialized to reflect a 
"blank" or cleared value. For each element in the 
temporary array, the element x, y coordinates are converted 
to distance or range from the radar site; if this range is 
within the radar data coverage area defined by circle 56, 
then a marker accompanying that temporary array element in 
the database is given an "unset" indication. For each 
radar data bin location of the radar product then being 
mapped, the bin x, y coordinate is converted to latitude 
and longitude; the resulting bin latitude and longitude is 
converted to a grid x, y location; the resulting bin x, y 
and grid x, y pair so found is assigned to an appropriate 
location in the look-up table for that site; the temporary 
array element x, y location is computed for that pair; and 
the temporary array element at that x, y location so found 
is given a "set" indication, showing that no hole is 
present at that location. Thereafter, the following steps 
are performed for each element in the temporary array 54. 
It is first determined whether that temporary array element 
is unset and thus is in the radar data coverage area but 
without a radar data bin assigned to it. If so, then the 
temporary array element x, y coordinates are converted to 
grid x, y coordinates; the resulting grid x, y coordinates 
are converted to latitude and longitud ; th resulting 
latitud and longitude valu s are c nv rted to radar bin x, 
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y coordinat s; and the resulting bin x, y coordinate and 
grid x, y coordinate pair are assigned to the look-up table 
for that single site radar product, to fill or to partially 
fill a hole in the mosaic grid. Finally, the resulting 
look-up table is assigned to a file in the database 46. 

The process of Fig. 9 to build a radar data mosaic 
using the look-up tables can be summarized as follows. 
First, for each location or element in the large grid 22, 
that grid element is initialized to a suitable 
predetermined value. Then, the following steps are 
performed for each radar site. The latest weather radar 
data product for that site is read in from the database. 
That radar product is then decoded, and individual data 
values from that product are assigned to a two-dimensional 
array. The corresponding prebuilt look-up table is then 
read in from the database 46. Then, for each entry in that 
look-up table, extract the corresponding bin x, y 
coordinate and its corresponding grid x, y coordinate from 
the look-up table. The product data value at that bin x, 
y location is then extracted from the decoded radar 
product. That product data value is then assigned to the 
corresponding grid x, y location obtained from the look-up 
table. These last three steps are repeated for each entry 
in that look-up table. Finally, the resulting combined 
radar data in the grid is encoded and saved to file in the 
database 46. 

Figs. 8A, 8B, 8C, 8D and 8E together show a structure 
diagram or flowchart for a method of building look-up 
tables. These look-up tables are essential components for 
the method of Fig. 4 and of its corresponding system for 
building mosaics of radar data in real time such as on a 
general purpose digital computer. The method of Figs. 8A- 
8E can be implemented as a computer software program that 
preferably can be run offline when time is not critical, so 
that time-consuming computations are performed in the look- 
up table building process and thus w 11 before th results 
of those computations would b needed. 
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As s ?wn in Fig. 8A, th program uses a manually 
created configuration file 101 that contains information on 
one or more sites for which a look-up table is to be built. 
Configuration file 101 contains the number of sites, and 
the latitude/ longitude and a site l.D. number of each site. 
At step 102, the program reads the number of sites from 
file 101 , and at step 103 reads in the site information 
from file 101 and stores that in the computer's memory for 
later use. Then, at step 104, the program enters a loop to 
process each such site on that list. This loop is executed 
until all sites are processed (determined at step 105), at 
which time the program terminates (step 106) . The 
processing by that loop for each site begins at step 107 
with the initialization of the temporary array by assigning 
the value of -1 to all elements of that array. 

Next, the temporary array is set up by the program, as 
shown by Fig. 8B, to indicate which elements of that 
temporary array are located within the radar coverage area. 
At step 110, this program loops for each x coordinate in 
the temporary array until all x coordinates are processed 
(determined at step ill). In this second loop, the 
temporary array x coordinate then being processed is 
converted at step 112 to the corresponding grid array x 
coordinate. At step 113, it is then determined whether the 
resulting grid x coordinate is outside of the boundary of 
the grid; if so, then the program skips to the next 
temporary array x coordinate (steps 110, 111 and 112). 
This could occur for sites located near the east or west 
edges of the grid (for example of the region of interest) 
where the temporary array 54 extends beyond the edge of the 
grid 22. Otherwise, for a valid grid x coordinate, the 
program proceeds to another loop starting at step 114 for 
each y coordinate in the temporary array until all y 
coordinates are processed (step 115). At step 116, the 
temporary array y coordinate then being considered by this 
third 1 op is converted t the corr sponding grid array y 
coordinat . If at st p 117 the resulting grid y co rdinat 
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is found to be outside of the boundary of the grid, then 
the program skips to the next temporary array y coordinate 
(steps 114, 115 and 116). That could occur for sites 
located near the north or south edges of the grid (or of 
the geographic region of interest) where the temporary 
array extends beyond the edge of the grid, otherwise, for 
a valid grid y coordinate, at step 118 the grid x, y 
position found by steps 112 and 116 is converted to 
latitude and longitude by using the projection equations 
that define the particular grid being used. For example, 
Lambert conformal projection equations could be used for 
this purpose. Since a grid is treated as being flat, but 
is used to represent a curved surface, some projection will 
be used to transform the grid into the curved surface or 
vice versa; various ways to accomplish such projection are 
well known to cartography. At step 119, the latitude and 
longitude is then converted to x and y coordinate values 
for the radar coordinate system, i.e. the coordinate system 
for the single site radar product. For step 119, an 
appropriate cartographic (or the like) conversion from 
latitude and longitude to the coordinate system of th 
radar product array is needed. The opposite such 
conversion was accomplished at step 118, except that as 
discussed above, the grid coordinates would not necessarily 
correspond to those of the single site radar product array. 
The x coordinate and y coordinate found at step 119 is then 
used at step 120 to compute the distance (in the radar 
product array) of the location defined by the x, y 
coordinate from the radar site. At step 121, it is 
determined whether the distance found at step 120 is within 
the radar coverage area such as is defined by circle 16 of 
Fig. 1; if so, then the value of the temporary array 
element, defined by the x, y coordinate utilized at steps 
112 and 116, is set to zero at step 122, to indicate that 
the corr sponding radar bin value is defin d by the singl 
site radar product array. Either way, the program g es on 
t the next y co rdinat at step 114. 
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Once th temporary array is set up by the m thod of 
Fig. 8B to identify the elements that fall within the radar 
coverage area, generation of the look-up table proceeds, as 
shown by Fig. 8C. As shown in Fig. 8C, generation of the 
look-up table begins at step 130, where the number of look- 
up table entries is set to zero. Starting at step 13 1, the 
program then performs a loop for each y coordinate in the 
radar product raster array until all y coordinates in that 
array are processed (step 132). For each iteration of the 
loop beginning at step 13 1, the program performs a loop 
starting at step 133 for each x coordinate in the radar 
product raster array; this loop is exited at step 134 to 
return to step 131 when all x coordinates in the radar 
product raster array are processed by the loop beginning at 
step 133. In the loop beginning at step 133, using the 
current product raster array x, y coordinate, the distance 
of that coordinate (in that array) from the radar site is 
determined at step 135. At step 13 6, it is determined 
whether the distance computed at step 135 is outside of the 
radar coverage area; if so, then the program skips to the 
next x, y position by returning to step 133. Otherwise, 
the position is within the radar coverage area and further 
processing is performed, as shown in Fig. 8C. In this 
further processing, at step 137 the x, y coordinate then 
being processed is first converted to latitude/ longitude. 
As used herein, coordinate or position refers to an x-value 
and a y- value that together (such as Cartesian coordinates) 
define a location in a system, array, grid, surface, or the 
like. At step 138, that resulting latitude/ longitude is 
then converted to the corresponding grid coordinate system 
x, y position. At step 139, the resulting grid x, y 
coordinate is checked for validity. If at step 139 it is 
found that the grid x, y position found at step 138 is 
invalid (not within the preset grid boundaries), then those 
coordinates do not belong in the look-up table then being 
constructed, and the program skips, to the next element in 
the radar product array b ing pr cessed, by returning to 
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step 133. Otherwise, at st p 140 the radar product rast r 
array x, y coordinate, and the corresponding grid x, y 
coordinate found at step 138 for that radar x, y coordinate 
are, assigned to the look-up table array then being formed 
in the computer's memory. At step 141, the total number of 
look-up table entries for that look-up table is then 
incremented by one; recall that before the program entered 
the step 131 loop, that number had been initialized to zero 
at step 130. After step 141, the grid x, y coordinat 
found at step 138 and entered in the look-up table at step 
14 0 is used at step 142 to compute the temporary array x f 
y coordinate corresponding to that grid x, y coordinate. 
At step 143, the value of the temporary array element at 
that x, y position is set to the value of plus one, to 
indicate that the corresponding grid box received at least 
one data value from the radar product raster array. Th 
program then returns to step 133. If at step 132 it is 
found that all of the y coordinates in the radar product 
array have been processed by the loop beginning at step 
131, then the program proceeds to the process shown at Pig. 
8D. 

After all elements in the radar product raster array 
have been processed by the method of, Fig. 8C, then th 
program locates all grid boxes within the radar coverage 
area that were not matched up with a coordinate or position 
from the radar product raster array, as shown in Fig. 80. 
This is accomplished by searching the temporary array for 
all unset elements given the value of zero at step 122 of 
Fig. 8B. In Fig. 8D, a loop beginning at step 150 is 
performed for each x coordinate in the temporary array 
until there are no more x coordinates to process (step 
151) . Within that loop, another loop beginning at step 152 
is performed for each y coordinate for the x coordinate 
then being considered by the step 150 loop until all such 
y coordinates are processed (step 153). In this other 
loop, aft r it is determined at step 153 whether there 
still remains any y coordinates t so proc ss, at step 154 
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the valu at th curr nt x, y locati n in a temporary array 
is evaluated. At step 154, it is determined whether the 
temporary array element at the x, y location then being 
addressed by the two Fig. 8D loops is equal to zero. If 
the value of that location is non-zero, then it is either 
outside of the radar coverage area or else it has been 
matched with a radar data bin location. Otherwise, when 
the value at that x, y location is still unset (in this 
example, assigned a value of zero) , then a hole exists in 
the mosaic. For each such hole, a match is needed to the 
radar data bin coordinate closest in position to that. If 
such a hole exists, then the grid x, y coordinate is 
determined at step 155 from the temporary array x, y 
coordinate considered at step 154. At step 156, the 
resulting grid x, y coordinate is converted to 
latitude/ longitude using the equations (discussed above) 
that define the grid projection. That latitude/ longitude 
found in step 156 is then converted at step 157 to the 
closest x, y coordinate in the radar product array 
coordinate system. At step 158, it is determined whether 
the radar product x, y coordinate found at step 157 is 
valid; if not, then the process returns to step 153. At 
step 159, it is then determined whether the grid x, y 
coordinate determined at step 155 is valid; if not, then 
the process returns to step 153. If the radar product x, 
y coordinate is valid and the grid x, y coordinate is valid 
(i.e., they fall within the region of interest), then at 
step 160 those coordinates are assigned to the lookup table 
array in the next available location in that array. Should 
an invalid such coordinate occur, the process then skips to 
the next temporary array element at step 153. After the 
radar product x, y coordinate and the corresponding grid x, 
y coordinate are added to the lookup table array, at step 
161 the total number of entries in the lookup table is 
increment d by on . From step 161, the process returns to 
st p 153. If n furth rye ordinate (s) is (ar ) pres nt, 
then the process returns t step 150 so that the n xt x 
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coordinate can be addressed in the manner described above 
for Fig. 80. If no more such x coordinates remain to be 
processed, then the program goes from step 151 to step 170 
of Fig. 8E. 

Once the lookup table is completed in the manner shown 
in Fig. 80, the program saves the table so determined in a 
file such as on the computer's disk in the manner shown in 
Fig. 8E. First, a unique file name for the lookup table is 
created at step 170 associating the current site 10 number 
with the file. At step 170, the file is then opened for 
output. At step 172, the program then builds a lookup 
table header in memory containing, among other information, 
the number of entries in the table, determined at steps 
130, 141 and 161. The program then writes to the new 
lookup table file 173 at steps 174 and 175. First, th 
header built at step 172 is written at step 174 to th 
lookup table file 173, followed at step 175 by the complete 
lookup table. At step 176, the lookup table file 173 is 
then closed. The program then returns to step 104 of Fig. 
8A to continue processing until lookup tables for all sites 
in the site list are complete. It is only necessary to 
rerun the program of Figs. 8A-8E if one or more sites are 
moved or added to the mosaic; real time running of the 
program of Figs. 8A-8E is not needed. 

Construction of a look-up table according to the 

present invention can be implemented in software. An 

example of such software, written in the Microsoft version 

of the C language, is given in the following pages. 

However, other programming languages can be used for this 

purpose. Figs. 8A-8E together generally show a flowchart 

for this software. 
/*«•*****«««*«*•*•***..«*•*.• * * * 

** UNISYS WIS MOSAIC PROCESSING SYSTEM 
** copyright 1994 Unisys Corporation 
#* 

** PROGRAM NAME: CR_LUT.EXE 

** DBSCRIPTI NtRead In th configuration file and build a lut for each eite. 
•» Maps Composite Rafl ctivity Product 36 and 38 to a 4x4ka National grid. 

*• DATA FILES NEEDED t 

*« C0NPI0.DAT ( Binary data file of sitae for which lut filee are 
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*+ to be built* Created by the conflg.exe program. 

** Located in the local directory.) 

** 

•* OUTPUT FILES CREATED i 

** CR4 ???.LUT (Maps NEXRAD Composite Ref product* to 4x4km grid. 

•* - There must be one for each NEXRAD site.) The ??? 

•* represents the NEXRAD site id number. 



/include <math.h> 
/include <atdio.h> 
/include <stdlib.h> 
/include <strlng.h> 
/include <limits.h> 
/include <memory.h> 
/include <dos.h> 
/include <fcntl.h> 
/include <lo.h> 
/include < s}*j\ types . h> 
/include <sys\stat.h> 

/define MAXWRITE 65534 

/define NODATA 0 
/define KIN NUM 0 
/define KIN~VALUE 0 
/define NO VALUE 0 
/define MIN X 0 
/define MINJf 0 

/define GRID X SIMjMCM 1280 
/define GRIDjCs IZB - 4KM 896 

/define LUT_MAX 190000 

/define GRID X HALF 4 KM 640 
/define GRIDjQ !AL 0 KM 448 

/define PR BINS_CR4 232 
/d fine PR~ROWS CR4 232 
/define TMP X ST«_CR4 250 
/define TMP~V~SI«_CR4 250 
/define THP_HALF_CR4 125 

/define BIN X_IDX 0 
/define BIN - T IDX 1 
/define GRID X IDX 2 
/define GRID~Y~IDX 3 
/define LUTJIiBth 4 

/define N PI 3.14159265356979323646 
/define M~PI 2 1.57079632679489661923 
/define HJIJ 0.78539816339744830962 

long lutjwy - {123456} ; 
typedef struct { 

long keyword; 

long num entries? 

short header sire; 

short site iH; 

short version; 

short date_ built; 

unsigned long time_built; 

short grid_res; 

short grid~wldthf 

short grid_height> 
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short bin_res; 
short bin~"width; 
short bin~height; 
short projection; 
double grid center_phi; 
double grid~center lands; 
double radar_phl; 
double radar_lamda; 

> lut_header_for£at; 

lut_header_format lut_header; 

typedef struct { 
double dx; 
double dy; 
double phi org; 
double laa3a_org; 
double phi value; 
double lamSa value; 
short *; 
short y; 
short pr id; 
short grld_id; 

> coords ; 

#def ine MAX_SITBS 200 

long site_xey - {24680} ; 
short nexrad_slte ■ {1>; 
short rrwds site » <2>; 
short yes "~{1>; 
short no • {0>; 

typedef struct { 
long keyword; 
short nuu_sitee; 

> site_heade reformat; 

site_header_format aite_f ile_header; 

typedef struct { 

short site id; 
short site~flag; 
short sita~build_lut; 
short site_use; " 
double site lat; 
double slte~lon; 
char aite_cKar [8 J ; 

> site_content s_f ormat ; 

aite_contents_f ormat aite_f ile_contents(KAX_SITBSJ ; 

double max_range » {464 .0}; 

double phi center_deg » 38.0; 
double iaxn3a_center_deg - -98,0; 

/• 

** Prototypes. 
•/ 

short rnd off (double fl); 
int 1 c km21atlon< coords *1 coords ); 
int 1~ c~latl n2km< coords *l"coords ); 
int r"c~lati n2grid( coords •! coords ); 
int 1 c~grid21atl n< coords *1 coords ); 
int nex~km21atlon< coords *l_coorde ); 
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int nex_latlon2km< coords *1 coords ); 
int nax latlon2bin( coords *T_coorda ); 
int nex~bin21atLon{ coords *l~coorde ); 
long huge_fwrite< void huge~*buf f er, long length, 
siae_t count, int fKw) ; 

/• 

** Main program. 
V 

int main() 
{ 

int status; 

unsigned int dos status; 

long nread; 

si*e_t nwrite; 

long Biit_writ«/ 

short binjc; 

short bin~y; 

short nuxD^eiftes; 

short j i 

long lut_size; 

char lut nana [100]; 

static sKort _huge lut [ LUT_MAX ) [ LUT_WIDTH J ; 
long entry; 

eh rt tmp_init - <-l>; 
ah rt tmp_unset » {0>; 
sh rt tmp sat ■ {1>; 
short goo3_status - {0>; 
double range; 
short grid_x; 
sh rt grid~y; 

d uble phi Hegi /* "Object* latitude (degrees). */ 
double lam3a_deg; /* -Object* longitude (degrees). •/ 
double phi radar deg; /* Radar latitude (degrees). */ 
double lam3a radar deg;/« Radar longitude (degrees). •/ 
static short" _huge tmp [ THP_X_S I ZK_CR4 ] [ TKP_Y_S 1 ZB_CR4 J ; 
short tx,ty; ~ 
double nx,ny; 

short gr id_radar_x , gr id_r adar_y ; 
c ords dummy coords; 
coords * locaT_coords ; 

FILB * fpr; 
int fhw; 



*• Initialise the local_coorda -pointer- so that 
** it points somewhere. - " 

*/ 

local_ coords - 6dummy_coords ; 



** Open the configuration file for "read binary*. 
*/ 

f pr-f open ( * conf ig . dat " , n rb" ) ; 
if (fpr— NULL) 

1 ( void ) pr int f ("Error opening the configuration file.\n"); 
return(l) ; 

y 



« Read in the configuration fil h ader. 
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nread-fread(6site file header, siteof < site header f ormat ) , 1 , f pr ) ; 
i f { nread-- NODATA ) ~ 

( 

( void) print f ("Error reading the configuration file header. \n")/ 
return(l); 

> 

/• 

** Check the configuration file header for the correct keyword 
** and exit the program if the keyword does not agree. 



if{ aite file header . keyword i- site key ) 
{ 

(void)printf ("Invalid configuration file keyword\n"); 
return ( 1 ) ; 

> 

/• 

** Extract the number of eitee from the configuration file header 
•/ 

numjsites * eite_f ile_header.num_sites; 
/* 

** Print the maximum number of sites. 
*/ 

(void) print f ("\nbuilding mosaic look up tables for %d eites\n\n" ,num_sites) ; 
/• 

** Read in the contents of the configuration file. The amount to be read in 
** is determined by the number of sites and the sits of the structure for 
*+ each site. 
*/ 

nread*fread(sita_f ile contents, num sitee*sizeof (site contents format ), 1, fpr) ; 
if ( nread—NODATA) 

{ 

(void)printf ("config.dat: fread (contents) error\n*)j 
return(l) ; 

) 

/* 

** Close the configuration file. 
•/ 

status - fc lose (fpr); 
if( status ) 

( void )printf ("Error closing the configuration file.\n"); 
return(l) j 

> 



** Loop through all of the sites and build a lut for each. 



for (j-MINJIUM; j<num_sitee; j++) 

if( (site file contents [jj. site flag « nexrad site) 66 
<site~flle~contente( j ) .site~buiid lut — yes) ) 

( 

/• 

•« Initial!* the temporary array with the following. 

The array is initialised by setting every array element. 
This array is used to determine where holes xist in the 
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*• initial mapping of NEXKAD data bins to the national mosaic. 

V 

for (tx-MIN X; tx<TKP X SIZE CR4 ; tx++) 
{ " - - - 

for(ty-HIN Y; ty<TMP Y SIZE CR4; ty*+) 
tmp[tx)tty]»tmp init; 

> 



/• 

** Extract the radar site latitude and longitude and assign the 

*+ values to simple variables and assign to the local_coords structure. 

phi radar deg- site file contents( j ] , site latj 
lamaa_radar_deg«site_f iTe_contenta[ j ] . site_lon; 

loca?l_coords->phi value ■ phi radar deg; 
local~coorde->lam3a_value - lamda_radar_deg; 

/• 

** Print the site id and radar latitude/ longitude. 
*/ 

(void ) print f("file# %3.3d ffSXKAD sitet %3.3d % j,eits file contents ( j J . site id); 
(void)printf <" radar lat/lont %f/%f (deg)\n", 
phi_radar_dag, lamda_radar_deg) ; 

/• 

»• Convert the radar latitude/ longitude of the current site to an x,y 

** coordinate in the national grid with the following routine 

** (which uses the Lambert con formal projection). The x and y 

*+ coordinates are stored in the simple variables* grid radar x, grid radar y. 

V " 

status*! c latlon2grld( local coords ); 
if (status)" 

{ 

( void) print f( "Error from l_c_latlon2gridr atata««%d\n" , atatus) ; 
return(l); 

> 

grid radar x • locel_coords->x; 
grld~radar^y - local~coords->y; 

/• 

** Next, unset the array elements which fall inside of the 

** range of the radar. This is done because this is the region 

•+ which must be checked for "holes" left by mapping NEXRAD bine to 

** the grid. First, loop through the temporary array in the x direction* 

*/ 

for (tx-MIW X; tx<TMP X SIZE CR4; tx++) 
/* 

** Convert the temporary array x coordinate to the corresponding 

** grid box x coordinate. 

•/ 

grid_x-(grld_radar_x - TMP_HALF_CR4 ♦ tx)> 



** Check that the x grid coordinate falls in the array. 
** Otherwis , the positi n is just skipped. 
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:; -oop throu , h all ef th . y bins in the teopo ^ ^ 

Jor(ty.MIN_Y; ty<THP_r_SI 2B _CR4; ty+> 
/• 

to th . 

grid_y. (grid _ radar _ y . THp _ HALp _ CR4 + ty); 

;; KM? £ loV^V^t^ «- -ay. 
if ( (grid_y>-MIN_y,« (orld _ y<ORID _ y _ SIZB _ 4iQ<) j 

statu. - l_c_grid21.tlon< local_coord. 

?l!ijl d *? " loc «l_coord.->phi valu., 
Wa^ - locIl_co 0 rd.r>Ii^r; alu<l; 

« 2 « < ET2.i& 1 2 1 £!: J"* 1 "?- portion 

loc.l_coord.->l aiB a a _ org - lS*«_rSdZr\tog, 
.t*tu. - n.x_latlon2km( local_coord. ^ 
/* 

*• t«nporary array valuVf? thf 5^? *" m l Vnmmt th « 
*• rang, of th. product di.tanc. i. within th. 

nx - local coords->dx; 
ny - local~coord.->dy; 

rang.-.qrt{ (nx«nx)*(ny*ny) ) , 
i £ ( r ang«<iiiax_rang. ) 

tmp[txj[ty)-trap_ un ..t/ 
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*« The look up table is built here. The first step is to "map" 

each NBXRAD bin to a grid x,y coordinate. The overall approach 

** la to convert the bin location to latitude /longitude and then 

** convert that latitude/longitude to a grid x,y coordinate 

** (which ie in the Lambert confonnal projection). The look up 

** table coneiete of a liat of entries, each containing a NEXAAD 

*• bin and row coordinate and the corresponding grid x and y coordinate. 

*• Later on, in the second step, holes in the grid are found and the 

** closest NSXRAD bin is mapped to the grid x,y coordinate of the hole. 
•/ 



** Initially, there are no entries in the look up table. 
*/ 

entry»NO_VALOK; 



*• Aajsign the radar latitude longitude value of the current aite 
** to the local coords structure. 

•/ 

local coords->phi org - phi_radar_degi 
local~coorda->lam3a_org - la»da_r adar_deg ; 



** Loop for each NBXRAD bin in the product. The outer loop 

** handles each row of data. 

•/ 

for(bin y-MIH V; bin y<PR ROWS CR4> bin_y++) 
< 

/* 

** The inner loop handles each bin in a row. 
•/ 

for(bin_x-MIN X; bin x<PR BINS CR4 1 bin_x++) 

< 

/* 

♦ • For the current bin, which ie in the MBXRAD 

•* coordinate system, determine the latitude longitude. 

•/ 

local_coords->x - bin_xj 
local~coords->y * bin_yj 

status-nex_bin21atlon( local_coorde )/ 

phi deg » local coords->phi_value; 
lamcTa_deg - locIl_coorde->lamda_vaiuei 

if (status—good status) 

/* 

** Por a valid x and y bin, take that latitude/ longitude 
•* and compute the corresponding grid x,y coordinate. 
•/ 

local coords->phi value » phi_deg; 
local"coords->iamaa_value - lamda_deg; 

statue»l_c_latlon2grid( local_coords )/ 

grid x - local_coorda->xj 
grid^y ■ local~coords->y; 
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if (statua»-good status) 

{ 

/* 

** For a valid latitude/ long itude, the corresponding 
** bin x,y and grid x,y coordinates are assigned to the 
*» look up table array and the number of entries 
** is incremented. 

V 

lut (entry) ( BIN X ZDX) • bin x; 
lut ( entry )(BXlTY~XDX] - bin~y; 
lut { entry ) { GRID X IDX) * grid x; 
lut(entry){CRID~Y~IDX) - grid^y; 

entry**/ 
/• 

** Then the grid x,y coordinate is used to compute the 
** temporary array x,y coordinate. If the temporary 
** array x,y coordinate is valid for the array, then 
** the array element is set. Any unset array elements 
** represent " holes* which will be handled later. 
•/ 

tx-grid x - grid radar x ♦ TMP HALF CR4; 
ty»grid~y - grid~radar~y ♦ TKP~HALF~CR4 ; 

if((tx>-MIN X)M<tX<THP X SIM CR4)K 

~ (ty>«MIN yyfc£(ty<TKP Y 5I«_CR4)) 

< 

tmp[tx] fty)»tmp seti 

> 

> 

} 

> 

> 



*• In the next step, go through the temporary array and find any unset 

** array elements which represent grid boxes within the product 

** coverage area that did not already receive a bin. For each of those 

grid coordinates, find the closest HEXRAD data bin and -map* it to 
«* the grid. This is where the holes are filled. Additional entries 
*• are made to the look up table for each set of KSXRAO bin, row and 
** corresponding grid x,y coordinates. 
•/ 

for <tx«KINJC; tx<THP_X_SIZB_CR4 ; tx*+) 

for<tyMIN Y; ty<TMP Y_SIZS CR4; ty++) 
{ 

/• 

** Check the current temporary array element to determine if 

** there is a hole with the following. (Otherwise, the loop just 

** goes on to the next element.) 

•/ 

if( tmp(tx}[tyl~tmp_unset > 

i 

/• 

** Convert the temporary array x,y coordinate to a grid x,y 
** coordinate. 

V 

grid x»(grid radar x - TMP HALF CR4 ♦ tx); 
grid^y- ( or id~radar~y - TMP~HALF_CR4 ♦ ty ) / 
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?atrtS.^: n | i r t i J d i: y COOrdiMt - *» the corr..po„di„g 

local_coord«->x ■ grid x; 
local_coord«->y ■ grid~y; 

statu. - l_c_grid21atlon< local_coord« ,, 

lo cal coorde->phi_value; 
lamda_deg - local_coord.->lamda_valu./ 

if(«tatu. «» good_atatua) 
/• 

11 f°f«? -*/f d 9rid *' v coordinate, convert the 

: : ^KT^tTO .<"'•••» 

local coords->phi_value - phi dee; 
local coorda->ian»da value - iindi d.c 
local coord 8 ->phi org - phi radar - d.<j; 
loc.l_cc*rd.->la»da_org - lamd."ad.1\i.g, 

etatue - nex_latlon2bin< local_coorde ); 

°ln_x - local coord.->x; 
bin_y - local2coorde->yj 

if (statu. « good_etatue> 

if((bin_x»MZN_X)ee(bin x<PR BINS CR4)M 

(bin y>-MIW r)**Tbin 5<pr rows CR4)cc 
griH x>-MlHjc)M (9r Ia x<orio J Bill 4HM)«e 
{ (9rid_y>-MIH_Y j **(grid>eORio;Y-SI8B-4IM!r 

/• 

11 U* 1 *^?"* bin coordinate, and the grid x,y 
coordinate, are valid, then aa.ign tho.e 

I. fn 0rdln * t- P * ir " to th *look up ?abl. "ray 
and increment the number of entrie.. Y 

lut[ entry JfBINJC IDX)-bin xt 
lut(entryJ[BIH Y~IDXJ-bin~y; 
lute entry) [ORIS 3f IDX)«grId x; 
lut ( entry J fCRID~Y _ iDX ) -grid~y ; 
entry**/ — *' 

/* 

° h ??£ } hm numb « r °* «ntrie. to determine if the 
oSt or ihm lo2? Print * and br "* 

•/ 

if (entry L0T_MJU) 



> 



S^l&jSllT™ ^ ' Ul1 ' •"^dNn-.entry,, 
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> 

lut_full:/ /* The goto branches her* to break from loop. Do nothing. */ ; 
/• 

** Create the look up table (lut) file name in which the site id number 

•* is embedded with the following. 

*/ 

{ void) sprint f ( lut_name, "cr4_%3 . 3d. lut\0" , site_f ile_contente[ j ) . site_id) ; 
/* 

** Print the lut file name. 
•/ 

(void) print f ( " look up table written to %s\n\n" , lut_name) ; 

** Create and Open the lut file for -write". 
•/ 

_fmode « 0_BINARV; 

f hw-creat ( lut_naioe , S I WRITE ) / 
if (fhw-— 1) 

{ void) print f( "Error creating lut file* %s\n" , lut_name) ; 
return ( 1 ) ; 

> 

lut header .keyword - lut_key; 
lut_header . num entries -""entry; 

lut_header.hea3er else - eizeof(lut header format); 

lut~header.site L3 - site_f ile_cont*nts( j l7aite_id; 

lut~"header .version • lj ~* 

lut~header.date built - 0; 

lut~header.time~built - 0; 

lut header. grid_rsa • 4; 

lut~header.gr id~width - GRID X SIZE 4KK; 

lut~header.grid_height - GRl5_7_SIZE_4KM; 

lut header.bin res - 4; 

lut~header.bin~width - PR_BINS_CR4i 

lut~header.bln~height - PR_ROWS_CR4 i 

lut~header .projection - li 

lut~header.gr id center_phi « phi_center_degi 
lut~header.gr id~center_lamda « lamda_center_deg; 
lut header. radar^phi - phi radar deg7 
lut~header . radar lamda - lamda_radar_degi 



** Write the lut file header structure to the lut file 

** with the following. 
*/ 

dos status - dos_write(fhw,fclut_header, <si*e_t)sizeof <lut_heade reformat) ,fcnwrit 
i f ( 3os_s t atu s) ~ 

(void) print f {"Error writing lut file header: %s\n" , lut_name) > 
return(l) ; 



f (nwrite— NODATA) 

( vo id) print f ("Error writing lut file header, nwrit - %d\n" , nwr its) ; 
return ( 1 ) ; 
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*« Compute the size of the lut file contents based on the number 
** of entries and numb r of bytes per entry. 

lut_size « entry*LUT_WIDTH*sizeof (short); 
/• 

** Write the lut file contents to the lut file 

** with the following. 

•/ 

size_write -huge_fwrite< lut , lut_size, 1, fhw) ; 

if (size writs-- NODATA) 
< 

(void)printf ("Error writing lut file contents. \n" ) ; 
return(l); 

> 

doe status-close ( fhw) ; 
if (3os status) 

{ 

< void) print f ("Error closing lut file* %s\n",lut name); 
return(l) ; 

> 

> 

> 

** Print a message indicating that processing is complete. 

(void)printf ("normal termination -> look up tables are complete \n\n") ; 

return (0) ; 

> 

/ *•»*.»**»*.•**»»*.•**»**.** 

* function: rnd_off( ) 

* Round off a real number and assign to a short. 
* •••••• * / 

short rnd off (double fl) 

< 

short 1; 

i- (fKO.O) ? (short)ceil(fl-O.S) : (short ) floor (f 1+0.5) ; 
return (i); 

> 

int nex km21atlon( coords *1 coords ) 

/• ~ 

* Names nex km21atlon 
* 

* Description i 

* Thie routine converts the x and y distance of an "object" in 

* the NEXRAO plane to the corresponding latitude and longitude 

* values in degrees. 
• 

•/ 

/• —————— .•—--input parameters - — */ 

/* 

* phi radar deg "Radar" latitude (degrees). 

* lam3a_radar_deg "Radar" longitude (degrees). 

* zx "Object" x distance from radar (in ka's) 
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where x increases toward* the east. 

zy "Object" y distance from radar (in Java) 

whsre y increases towards the north. 

/ 



# .„ — ——-output parameters— — — — — •/ 

* 

phi deg "Object" latitude (degrees). 

Iam3a_deg "Object" longitude (degrees). 

* , 



/* -Internal Variables- */ 

double phii /* "Object" latitude (radians). */ 

double lamda; /* "Object" longitude (radians). •/ 

double phi_radar; /* "Radar" latitude (radians). */ 

double lamda, radar; /* "Radar" longitude (radians). */ 

double radlul - {6380.}; /* Average earth radius (in kilometers). »/ 

double constant - {135.}/ /* Constant. */ 

double deg2radians ■ {MJPI/180.); /*Converte degrees to radians. */ 

double sp; /• Intermediate term. */ 

double *ain_theta; /* Intermediate term. */ 

double xcoe~theta; /* Intermediate term* */ 

double *sin~s; /* Intermediate term. */ 

double rcoe~s; /* Intermediate term. */ 

double seinjphi; /* Intermediate term. */ 

double zcoejphi; /• Intermediate term. */ 

double lamda delta; /* Intermediate term. */ 



double *x; 
double zy; 

double phi radar_deg; 
double lam3a_radar_deg; 
double phi deg; "~ 
double lam3a_deg; 

/* 

** Extract the radar latitude/longitude (in degrees) from the l_coords 
* * structure . 

*/ 

phi radar deg • 1 coords->phi_org; 
lam3a_radar_deg •~l_coorda->lamda_org; 



** Extract the object position (in kilometree) from the l_coorde 

** structure. 

•/ 

zx • l_coords->dx ; 
zy » l_coords->dy; 

/* Check the x and y distances from the radar location. */ 

if( ( (*x>-0.1)«(*x<0.1> > && ((zy>-0.1)M(*y<0.1)) ) 
< 

/• 

* This handles the case when the x and y distances are both 

* near sero (the object is essentially at the radar). 

* Set the bject latitude and longitude (in degrees) to the radar 

* latitude and longitude (in degrees). 
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phi deg - phi radar deg; 
^ lamda_deg - lamda_ridar_deg, 
alas 

{ 

/* 

* th. object latitudT and - th * C°*PUt. 

Cd£ a SSJ ^l£ dar -2' 9 * d *92radiana, 
xaawa_radar - lamda_radar_dag • dag2radiana, 

/• comput. th. diatanc. to th. obj.ct (in >cilom.t.r., . 

*P - egrt( (zx«»zx) + (zy*zy) ); 
/• Comput. th. .i„ and co. of th. angl. to th. obj.ct. 

z.in_th.ta - (zx/zp), 
zcoa_theta - (zy/zp); 

isssrjx jrsra: zsjiss ra-s^-saar- 

" An -" " !r.d?£; adAu "J/ ( 1 -°-<< «»■*** . «p ,/ 

(radiua * radius))) ); r f/ 

zco._a - .qrt( l.o - (zain_a • «. in , ; 
/• 

• SKS.?" int «»*<*^t. .in and oo.in. t.rm. ox th. obj.ct 

*/ 

z.in_pht - (sin(phi_radar)»zco. • ♦ 

co. (phi_radar) •z.in_.»acoa_th.ta) ; 

zcoa_phi - .qrt( 1.0 - (zain_phi*z.in_phi) ); 
/• Co«put. th. diff.r.nc. b.tw. n obj.ct and radar longitud.. ./ 
la»da_d.lta - ..i„ ( „ in _. . „ in _ tnata , aeoajphL 

/• comput. th. obj.ct latitud. and longitud. U „ radian.,. 

i»LZ atan2 i «-in_phl,zco. phi ,, 
lamda - lajada_radar+lamda_d.lta; 

/• convrt th. latitud. and longitud. from radian, to d.gr—. ./ 

Phi d.g - (phi/d.g2radian.); 
lamaa_d.g - (lamda/d.g2radian.) ; 

It«I n ur.T latitud -A°ngltud. dn d.gr—> to th. l_coord. 
coord.->phi_valu - phi_d.g; 
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l_coord»->lamda_value ■ lamda_d 9; 
I* 

Return a zero value. 

•/ 

raturn(O) ; 

} 

/ * • * — / 
int nex_latlon2km( coords *l_coords ) 

/• 

* Name: nex_latlon2km 
* 

* Description x 

* This routine converts the latitude and longitude of an "object" in 

* the NEXRAD plane to the x and y components of distance 

* in kilometers « 
•/ 

< 

/* — — — —Input parameters™——— — — */ 

double phi_deg; /* "Object" latitude (degrees). */ 

double lamda deg; /* "Object" longitude (degrees). */ 

double phi radar_deg; /* "Radar" latitude (degrees). */ 

double lamda_radar_deg; /* "Radar" longitude (degrees). »/ 

/* —-Output parameters— — — •/ 

double *x; /* "Object" x distance from radar (in km's) */ 

/* where x increases towards the east. */ 

double zy; /* "Object" y distance from radar (in km's) */ 

/+ where y increaeee towards the north. •/ 

/* Internal Variables */ 

double phi; /* "Object- latitude (radians). */ 

double lamda; /• "Object" longitude (radians). «/ 

double phi radar; /* "Radar" latitude (radians). */ 

double lamHa radar; /* "Radar" longitude (radians). */ 

double radius - <6380.>; /* Average earth radius (in kilometers). */ 

double constant « {135.}; /* Constant. */ 

double deg2radians - <M_PI/1B0.>; /'Converts degrees to radians. */ 
double zterma; /* Intermediate term A. •/ 
double ztermb; /* Intermediate term B. */ 

double zsin s; /* Intermediate termt sin of angular great circle. */ 

double zcos~s; /* Intermediate termt cos of angular great circle. */ 
double zd; ~ f* Intermediate term D. */ 

/• */ 

phi deg ■ l_coords->phi value; 
lam3a_deg - l_coorda->lamda_value; 

phi radar deg ■ 1 coorda->phi_org; 
lam3a_radar_deg «~~ l_coor de -> 1 amda_o r g ; 

/* Convert the radar latitude from degreee to radians. */ 
phi_radar - deg2 radians * phi_radar_deg; 

/• convert the radar longitude from degrees to radians. •/ 
lamda_radar - deg2radians * lamda_radar_d g; 

/• Convert the object latitude from degrees t radians. •/ 
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phi » deg2radians * phi_d g; 

/* Convert the object longitude from degrees to radians. */ 
lamda • deg2radians * lamda_deg; 

/* Compute intermediate terme a and b. */ 

zterma - coa(phi) * ain( lamda - lamda_radar ); 

ztermb ■ ( ( co« ( phi_r adar ) * ain(phi))- 

(sin(phi_radar) * coe(phi) * coa(lamda-lamda_radar ) ) ); 

/* Compute intermediate sine and cosine of the angular great circle 
distance from the radar position to the object position. */ 

zain_s - sqrt ( ( zterma* zterma) +( ztermb* ztermb) ) ; 

zcos_s - sqrt( 1 - (zsin_s*z«in_s) ); 

/• Compute intermediate term D. »/ 
zd « (constant • zsin_s) + radius; 

/* Compute the x distance (in kilometers). */ 
zx - zd * coe(phi) * sin(lamda-lamda_redar) ; 

/* Compute the y distance (in kilometers). */ 

zy - zd * (sin(phi) - ( ain <phi_r adar ) *zcoa_e) ) / cos (phi_r adar ) ; 

l_coords->dx • zx; 
l_coords->dy » zy; 

/• 

** Return a zero value. 

•/ 

return (0) ; 

> 

/ * • "* 

int nex_latlon2bin( coords *1 coords ) 

( 

short bin_x, bin_y; 
int status; 
double zx,zy; 

/* 

** Convert the lattitude/longitude to kilometers in the NBXRAO 

** coordinate system with the following. 

•/ 

status - nex latlon2km( 1 coords ); 
if (status) " 

{ 

return(l) ; 

> 

zx - l_coords->dx; 
zy - l~coorde->dy; 

/• 

Convert the x and y compon nts of distance from kilometers to 
*• date bin coordinates. Round off the floating point values as 
•* they are assigned to integer variables* 

•/ 
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bin x » rnd off( ( *x+462 . ) /4.0 ); 
bin^y m rnd~off( (462.-sy)/4.0 )i 

/• 

** Assign the x and y bin coordinates to ths 1 coords structure. 
•/ 

l_coords->x - bin_x; 
l~coords->y - bin^y/ 

/* 

•* Check that the x,y bin coordinate is within the valid range 
•* which determines the return value, A nonzero return value 
*• indlcatee that the latitude longitude corresponds to a 
** nonvalid x,y bin coordinate. 

V 

if ( (bin X>«MIN X)&&(bin x<PR_BINS CR4)M 

Jbin_y >-MIN_r ) && ( bin_y<PR_ROWS_CR4 ) ) 

/* 

*» This handles valid x,y bin values. A zero value is returned. 
•/ 

return(O) ; 

> 

return<2) ; 

} 

/ *** **** •*••/ 

int 1 c latlon2grid< coords *1 coords ) 

< 

short grid_X/ grid_y; 

int statue? 

double zx,zy; 

double boxaize - {4.0>; 

/* 

•* Convert the incoming latitude longitude to a distance in 

** (nominal) kilometers using the Lambert conformal projection. 

•/ 

status » 1 c latlon2km( 1 coords ); 
if (status) 

< 

/• 

** If a nonvalid status occurs then return from this routine 
** with a nonsero status. 

•/ 

return(l) ; 

> 



** Extract the x,y components of distance from the 1 coords structure. 

*/ 

«x * l_coorde->dx ; 
zy » l~coords->dy ; 



** Convert the x,y components of distance to grid x,y coordinates. 
** The valuee are rounded off as they are converted from floating 
** point to integer. 
•/ 

grid x - rnd f f ( zx/boxsize)+ORXD X HALF 4 KM; 
grid^y • GRl5_Y_HALF_4XM-rnd_of f (Iy7boxsIzs) / 
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/• 

** Assign the grid x,y coordinates to the 1 coords structure 
** with the following. 

l_coorde->x • grid x; 
l_coords->y ■ grid^y; 

/* 

** Check that tha grid x,y coordinata values rapraaant a valid 
** array index. Nonvalid values result in a nonzero return. 
•/ 

if ( (grid x>«HIN X)**<grid x<GRID X SIZS 4KM)6& 

(grId_y>-MTN_Y)66(grId_y<GRlo;y - SIZB - 4KM) > 

/• 

** This handles the valid grid x,y valuea. A zero value is returned. 



return(O) ; 

> 

return(2) / 

> 

/ •..«•*».«.., 

int nex_bin21atlon( coorda *l_coorde ) 

short bin x, bin_y; 
int statue; 
double zx,ty; 
doublet range; 

/* 

** Extract the NHXRAD data bin coordinates from the 1 coords structure. 
•/ 

bin_x - l_coorda->x; 
bin~y - l_coords->y; 

/* 

Check the x,y coordinate pair for a valid array index value. 
* * A non valid coordinate reulte in a non zero return from the routine. 
•/ 

i f ( ( binjc>«MIW_X ) fifi ( bin_x<PR_BINS_CR4 ) &fi ( bin_y>»MIN_Y ) « ( b in_y < PR_ROWS_CR4 ) ) 

/* Convert the bin x,y coordinate to the x,y componenta of diatance 

** where the reeult ia a pair of floating point valuea. 

V 

xx - ((bin x*4.0)-462.0) / 
«y " <-<<bin_y*4.0)-462.0)) j 

/• 

** Use the components of dietance to compute the range in kilometers. 

./ 

range - eqrt( <*x*zx)+( xy**y) ); 
/• 

*• Check the range to insure that it ia in the NBXRAD coverage area. 

** otherwls a nonzmr value ia returned. 

•/ 

lf(range<aax range) 
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*+ This handle* the coordinates located within the NEXRAD coverage 
** area. Assign the x,y components of distance to the l_coords array. 



l_coorde->dx » zx; 
l2coorda->dy • zy; 

/* 

** Compute the latitude/longitude of the point with the following 
*• routine which converts x,y distances in the NEXRAD plane to 
** latitude longitude (in degrees). 



status ■ nex km21atlon( 1 coords ); 
if (status) ~ 

< 

/* This handles a nonzero status returned indicating non valid 
** coordinates were passed. Return a nonzero status value. 
*/ 

return ( 3 ) ; 

> 

•* This indicates that valid latitude/ longitude values 
** were computed. A zero value is returned* 



return (0) / 

> 

return(l) ; 

> 

return(2) ; 

> 

/ .eeeeeeeee / 

int 1 c grid21atlon( coords *1 coords ) 

{ 

short grid_x, grid y; 

int status? 

double zx,zy; 

double boxsize * {4.0}; 

/• 

** Extract the grid x,y coordinate from the 1 coords structure. 
*/ 

grid^x « l_coorda->x; 
grid^y * l2coorde->y/ 

/• 

** Check that the x,y grid coordinate is a valid array index. 

** Otherwise, a nonzero status is returned. 

•/ 

if ( (grid x>-MIN X)fte(grid x<GRXD X SIZE 4KK)fi& 

(grid y>-MIN Y ) &6 < grId~y<ORID * SIZE 4 KM) ) 
{ ~ 
/* 

** This handles a valid x,y grid coordinate, convert the x,y grid 

coordinate to x r y distance components in (nominal) kilometers. 
*• The y ordinate is reversed. Then assign the distance components 
** t the 1 coords structure. 

V 

zx - (grid x * GRID X HALF 4RM) * boxsize; 
zy • < GRID jr_KALy_4E>T- grld_y) •boxsize i 
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l_co rda->dx - zx; 
l~coords->dy - xy; 

/• 

** Convert the x,y component s of distance (in kilometers) to 
** latitude/longitude using the Lambert conformal projection 
** with the following routine. 
*/ 

atatus«l c km21atlon( 1 coorde ); 
if ( status )~ 

i 

/* Thie handles the caee when a nonvalid condition has occurred. 
** K non tero value is returned. 

V 

return (2 ) ; 

> 

/* 

** This handles the valid condition and a zero value is returned. 
•/ 

return (0) ; 

} 

return(l) ; 

) 

/••••• •••••• ******** *** ********** 

* functions 1 c latlon2km( ) 

* Convert 1 at i£uae/ longitude to x/y using Lambert Conformal. 
********* *************** * *********** / 

/* Convert latitude /longitude to x/y •/ 

int 1 c_latlon2km( coords *1 coords ) 

i 

static double deg2radians - {M_PX/1B0.}; /* Converts degrees to radians. */ 

static double phi 1 deg * <33>; /• Standard latitude 1 (in degrees). •/ 
static double phi~2~deg « {45); /• Standard latitude 2 (in degrees). •/ 

static double radius * {6380. >; /* Average earth radius (in tan's). */ 

static double phi 0 deg - {38. 0); /* Grid origin latitude (in degrees). */ 
static double laia3a~0_deg ■ {-98. 0>; /* Grid origin longitude (in degrees). */ 

/* 

** The following constants were calculated for the grid center , earth radius, 
*• and standard latitudes defined above. DO NOT use this routine for other 
•* centers without recalculating the constants. 

•/ 

static double phi 1- {0.575958653158129}; /• Standard latitudes (in radians). */ 
static double phi~2- {0.785398163397448); /* Standard latitudes (in radians). •/ 

static double phi 0» {0.663225115757845}; /* Grid origin latitude (in radians). * 

static double lam3a_0- {-1.710422666954443}; /• Grid origin longitude (in radians). 

static double rho 0- {7931.823624772932817}; 

static double n- 70.630477697315427}; 

static d uble f» {1.955000201593793}; 

static double rf-{12472. 901286168398656); /• Intermediate term »(radius*f) */ 

double phi deg; 
doubl lam3a_deg; 
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double zx; 
doubl zyj 

double rho; 
double theta; 

double phi; /• current latitude value (radians). */ 

double lamda; /* Current longitude value (radians). */ 

/* 

** Extract the latitude and longitude of the object (in degrees) from 
** the 1 coords structure. 

phi deg » l_coords->phi value; 
lamda_deg -~l_coorda->lamda_value; 

/* 

*• Convert the latitude and longitude of the object from 
** degrees to radians. 

V 

phi - phi_deg*deg2radiana; 
lamda » lamda_deg»deg2radians; 

rho - < rf/(pow(tan(M_Pl_4 + (phi/2)), n)) ); 

theta - n*( lamda - lamdaj))/ 

/* 

** Compute the x and y components of distance (nominally in kilometers), 

zx - (rho^sin(theta)); 

zy - (rho_0 - rho*cos (theta) ); 

/• 

** Assign the x,y distance components (nominally in kilometers) 

** to the 1 coords structure. 

•/ 

l_coords->dx * zxj 
l2coords->dy * *y ; 

/• 

*• Return a sero value. 

•/ 

return (0) ; 

> 

lnt l_cJem21atlon( coords *l_coords ) 



ii 1 c km21atlon 



* Description* 

* This routine converts the x and y distance components of an "object" 

* in the Lambert Con formal projection to the latitude and longitude 
in degrees. 



* 

•/ 



/• *zx "Object- x di«tanc« P frcm radar (in km' a) •/ 

/* whara x incraaaaa towards tha aaat. •/ 

/• •/ 
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/* *ry "Object* y distance from radar (in km'e) */ 
/* where y increaaee toward* the north* */ 
/• / 

/* — — — — — -Output parameter a*-- — ~— — */ 

/* phi deg -Object" latitude (degrees). •/ 
/* lam3a deg "Object* longitude (degrees) . •/ 

/• — V 

{ 

static double phi 1- {0.575958653158129}; /* Standard latitudes (in radians). * 
static double phi~2« {0. 785398163397448); /* Standard latitudes (in radians). * 

static double phi 0» {0.663225115757845); /* Orld origin latitude (in radians), 

static double lanflaj)- {-1.710422666954443}; /* Grid origin longitude (in radians) 

static double rho_0« {7931.823624772932817}; 
static double n- {0.630477697315427); 
static double f- {1.955000201593793}; 

static double rf«{12472. 901286168398656); /* Intermediate term *(radius*f) */ 

static double radius » {6380.}; /* Average earth radius (in kilometers). */ 
static double deg2radlans ■ {K_PI/180.}; /* Converts degreee to radians. */ 

static double phi 1 deg - {33}; /* Standard latitude 1 (in degrees). */ 
static double phi~2~deg ■ {45}; /• Standard latitude 2 (in degrees). */ 

static double phi 0 deg - {38.0}; 
etatic double lam3a~0_deg ■ {-98.0}; 

/* _ Internal Variables——-——---- •/ 

double phi; /* "Object" latitude (radians). */ 

double lamda; /* "Object" longitude (radians). */ 

double rho; 
double theta; 

double ax, zy; 
double phi deg; 
double 1 am3a_deg ; 



/* < 

/* 

** Extract the object distance (in km's) from the l_coords structure. 
./ 

zx • l_coords->dx; 
zy » l~coords->dy; 

rho - sqrt( (zx*zx) * ( (rho_0-zy>* (rho_0-zy) ) ); 
/• 

«« compute the latitude in radians. 
•/ 

phi - {<2*atan( pow( (rf/rho) , (1/n) ))) - M_PI_2); 

theta * atan2( zx,(rhojO - zy) ); 

/• 

** Compute the longitude in radians. 
•/ 

lamda - ( (theta/n) ♦ lamdaj)); 
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•* Convert the Xatitud /longitude from radians to degrees. 

•/ 

ph i deg-ph i / deg2 r adiana ; 
lamda_deg-lamda/deg2radiane t 

/• 

** Assign the latitude and longitude (in degrees) to the 1 coords structure. 
*/ 

i_coords->phi value « phi deg; 
l~coords->lam3a_value - lamda_deg; 

/• 

** Return a xero value. 
•/ 

return (0) ; 

> 

** Function to write contents of huge arrays to a file. 



long huge_fvrite( void huge *buffer, long length, 
sise t count, int fKw) 

/• 

** Description: Function to write contents of huge arrays to a file. 

•/ 

< 

char huge *pnt_buffer; 
long remain; 
slM_t amount; 
sixe^t nwrite; 
long~total; 

unsigned int dos_status; 

pnt — buffer » buffer; 
remain- length ; 
total«0; 

while ( remain>0 ) 

{ 

if ( remain>MAXWRXTS ) 
< 

amount •MAXWRXTK ; 
remain — (long) amount; 

dos status - __dos write (fhw,pnt_buf far, amount ,finwrite) ; 
if (3os_status ) ~ 

* (void)printf (*dos_write bad status\n"); 
raturn(O); ~ 

> 

pnt_buffer amount; 
total nwrite; 

if (nwrite « 0) 

{ 

return (total); 

> 

> 

Is 

amount- <sire_t> remain; 
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remain-G ; 

doa atatua m _doa_write(fhw,pnt buf for, amount ,&nwrite) ; 
if(3os status) 

< 

(void)printf ("dot write bad ■tatu«\n")/ 
raturn(O) ; 

> 

pnt_buffer +• amount; 
total (long)nwrite; 

if(nvrite — 0) 

{ 

return (total ) j 

} 

y 

> 

return (total) j 
> 

Figs. 9A, 9B, ?C, 9D and 9E together show a structure 
diagram or flowchart depicting a program capable of 
building mosaics of radar data. The program of Figs. 9A-9E 
is designed to run in real time by using the lookup tables 
built by the offline program of Figs. 8A-8E. 

The program of Figs, 9A-9E, when started, runs as a 
continuous process. As shown in Fig. 9A, a manually 
created configuration file 200 containing a list of sites 
to be used in the mosaic is used to initialize the process 
at steps 201 and 202. First, the process reads in, from 
the configuration file 200, the number of sites (step 201) 
and the site list (step 202) that is then stored in memory 
for later use. Next, a loop is entered (step 203) that 
runs continually (shown in Figure 9B as a so-called "DO 
FOREVER" loop) . Each execution of the loop results in the 
creation of a new mosaic that is a compilation of the most 
current data from the individual sites specified in the 
site list. In preparation for the creation of the mosaic , 
all elements of the two-dimensional grid within the 
computer's memory are set (step 204) to the same 
predetermined value such as zero (i.e., the grid is cleared 
of any existing data). Next, an index "to the current site 
in the site list" is set (step 205) to point to the first 
site n th list. Then, another loop is entered that 
processes (step 206) the radar sit s on th sit list one 
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Figure 9D shows how, once the look-up table is read 
in, the single site radar data (which resides in the 
product raster array) is assigned to the two-dimensional 
grid. A loop is entered to process all of the entries in 
the look-up table (step 260). When all processing is 
complete, the program goes on to the next site (step 261). 
First, the x and y indices into the two-dimensional product 
raster array are extracted from the look-up table (step 
262) and checked for validity (step 263). (An invalid 
entry should never occur in practice, but if it should 
occur the program skips that entry and goes to the next 
entry.) The product data value is then extracted (step 
264) from the product data array. Then, the x and y indices 
into the two-dimensional mosaic grid are extracted from the 
look-up table (step 265) . Those coordinates are also 
checked for validity (step 266) . (Again, an invalid entry 
should never occur in practice, but if one should be 
detected then the program skips to the next entry.) Next, 
the product data value is checked for validity (step 267) 
and any invalid data values cause the program to skip to 
the next look-up table entry. Then, the product value is 
compared (step 268) with the existing grid value and, if 
greater, is assigned (step 269) to the grid at the x,y 
location, and, if not greater, is skipped. Processing 
continues in that way for all remaining entries in that 
look-up table. 

After all available sites have been so processed, the 
data in the grid is converted to a new product as shown by 
Figure 9E. First, the values in the grid are converted 
(step 280) to the output product format and assigned to a 
buffer in the computer's memory. Then, a header is built 
(step 281) , also in the computer's memory, containing the 
current date and time as well as the product data length. 
A unique output filename is created (step 282) and a file 
of that name is op ned for output (step 283). Next, th 
data is actually written to the output product fil 284. 
First, the product header is written (step 285) follow d by 
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the product data itself (step 286). Finally, the output 
product file 284 is closed (step 287) and the program goes 
on to start the next mosaic. 

Direct mapping of individual radar data values to a 
mosaic grid, using one or more prebuilt lookup tables 
constructed with the software given above, can be 
implemented in software. An example of such software, 
written in the ANSI version of the C language usable on a 
Sun Microsystems workstation is given in the following 
pages. However, other programming languages can be used 
for this purpose. Figs. 9A-9E together generally show a 
flowchart for this software. The software appearing in the 
following pages implements steps 48 and 50 of Fig. 4. The 
software shown above, and the software shown below, 
together implement steps 34, 36, 38 and 40 of Fig. 3 



** 
*» 
** 

*« 
*• 
** 

• * 
*« 
•• 
*« 
•• 

• * 
** 



UNISYS WIS MOSAIC PROCESSING SYSTEM 
Copyrigth 1994 Unlaya Corporation 

PROGRAM NAMEl BUILD. EXE 

DESCRIPTION! Raad in individual radar product* and corraaponding 
took Up Tablaa and map tha product data to a grid. 
Than taka tha raaulting grid valuaa and produca a 
naw product. 

ASSUMPTIONS t Products for tha daalrad aitaa raalda in %hm apaclfiad 
dlractoriaa. Look Up tablaa for tha daalrad aitaa axiat 
in tha "lut* diractory. 



#lnclude <math.h> 
#include <arrno.h> 
#includa <tixaa.h> 
lincluda <atdlo.h> 
#include <atdlib.h> 
^include <atring.h> 
lincluda <ctypa.h> 

/* 

*• Prototype of aubroutina. 

int gat_lataat < ahort , ahort , char *)* 

fdsflna BIN X IDE 0 
#daf ina BIN~Y~IDX 1 
#daf ina GRID x IDX 2 
#defina ORID~Y~IDX 3 
#d fine LUT_Wl5TH 4 

#daf ina MAX DATA SIZE 300000 
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*d*f ine NUM_DATA_LSVSLS 16 

#def ine NOD ATA 0 
#define MIN NUM 0 
/define MI UPVALUE 0 
#define NO VALUE 0 
#define MlRjC 0 
#define MINJf 0 

/define GRXDJC SZZS 1280 
/define QRIDJCJSXZZ 896 

#def ine lut max 50000 
/define PR BINS 232 
/define PRJtOWS 232 

long lut_key « {123456} ; 
typedsf struct { 

long keyword; 

long nua_entries; 

short hiidsr airs; 

short site i3; 

short version* 

short dat*_built; 

unsigned long tiM_bullt; 

short grid rss; 

short grid~width; 

short grld_hsight; 

short bln_rss; 

short bin_width; 

short bin~hsight; 

short projection; 

doubls grid centsr_phi; 

double grid~center_lamda; 

doubls radaxjphi; ~ 

doubls radar_lamda; 

> lut_header_format; 

lut_header_format lut_header; 

typedef struct { 
doubls dx; 
doubls dy; 
doubls phi org; 
doubls lam3a_org; 
double phi value; 
double lam3a_value; 
short x; 
short y; 
short pr id; 
short grld_id; 

> coords; 

long sitejiey - {24680); 
short nexrad site • {!}; 
short rrwds_slte - {2>; 
short yes -~<1) ; 
short no ■ <0); 

typedsf struct { 
long keyword; 
short nun sites; 

> eite_headar_format; 

slte_header_format sit _f ilsjisadar; 
typed** struct { 
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short sits_id; 
short aits flag; 
short sits3?uild_lut/ 
short sits_u»s; 
double sits_lati 
doubls sits Ion t 
} site_contenta_f ormat ; 

#defina MAX^SITES 200 

aita_contents_f ormat sits_f ila_contsnt»{MAX_SXTBS] ; 

PILK * fpwj 
FILB * fpr; 
FIL* * fprj?rod* 

char * buf; 
long data_eixa; 
long nwritsy 
long nroad; 

static Int errno - 0; 

typed* £ struct { 

short sits id; 

char *aite~textj 
> s it e _path_f ormat * 

#def ins HUM SITS PATH 30 



Example sits path information. 

•/ 

static sitsj>ath_format a it a_path { NUX_SITE_PATH ] « 

354, «AL maxwell AFB/", 
395, -AlTlittleroek/- , 

524 , -A3Q>hoenix/ " , 
347, -CO denver/*, 
351,*DB~dover/ w , 

302 , *rL~melbourna/ " , 
728,-FZTmlaAi/'', 
337,-2lTchicago/", 

381 , • IH~indianapol is/ • , 
350 , *KS~dodg*city/ " , 

366, "KS goodland/- , 
554,*KS**topeXa/*, 
562 , •KS~wlchlta/* , 
349,"MI~detroit/-, 
308 , •MO_st louls/ ■ , 
342, -MS columbus/-, 

382, *MS~jackson/-, 

367 , -NaThastings/ • , 
372,-mTgriffia AFB/" , 
340 , *OH"cievela5d/ ■ , 
305 , -OK~f rederick/ " , 
301 , •OKfnorman/", 
557,"OK""tulsa/", 

525 , -PA^pittsburgh/ * , 
313, *TX amariilo/-, 
332,*TX~cantral/", 
353 , *TX~dyess_AFB / - , 
378,*TX~houst n/", 

303 , •VA^stsrling/ • , 
000,-UOTIFIHKO/-, 
>l 
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typedef struct { 

short prod_id; 

char *procTt«xt; 
} prod_path_f ormat ; 

#def ine NUM_PROD_PATH 1 

static prod_pa t h_f ormat prodjpath(NUM PROD_PATH) 
{ 

38,"COKP REF/* , 

>l 

/• 

** Main routine to build mosaic. 
•/ 

roain( ) 

{ 

^define initial_x 0 
#define init*lal_y 0 
int statua; 

unsigned short * index; 
unsigned short *st_index; 
unsigned char bytej 
short n_bytee; 
short run, col; 
char *cp; 

unsigned short lengthl, length2; 
short k,m; 
short p; 
short bin; 

unsigned int timely time2; 

tlM_t tot sec; 

unsigned aEort opcode; 

short msg_code; 

short input_eource_ld; 

short output_sourcs_id; 

short op_ mode ; 

short vcp; 

short nrows; 

long lut_siae; 

long age~minutee ; 

static long max_age_ininutee • {15>; 

long age_ seconds; ~ 

short grld_x; 

short grid~y; 

short num_slteei 

short site_id; 

short prod~id; 

short j; 

short header size - (136); 

short product (PR BINS J [ PR_ROWS ] ; 

unsigned short g?idtOR2DJt_SIMHGRIDjr_SXK) ; 

short row; 

short box_x,box_y/ 

short row^index; 

short bin^lndex; 

short max~"col; 

short numEer_rows ; 

short number"~bins; 

short start Bin; 

short start~row; 

char lut na5e(S12); 

char pr©3 name (512); 

char prod~file(25J; 

char *s; ~~ 

sh rt lut(Ltrr_KAX) [LUTJTCDTH] ; 
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long entry; 

long num_antry; 

short preV_col; 

short nbytea; 

long layer length; 

long layer? length; 

long block_Iength; 

long meesage_length; 

long total length; 

short scale_x, sealery; 

unsigned short *nbytes_ index; 

unsigned short *mew_layer_len_index; 
unsigned short *lsw~layer~len~index; 

unsigned short *nisw_layar2_len_index; 
unsigned short *law~layer2_len~index; 

unsigned short *raswJblock_len_index; 

unsigned short *lsw~fclock~len~index; 

unsigned short *msv~msg_len_index ; 

unsigned short *lswjnsg__len~index; 

unsigned short snsw layer_len; 
unsigned short lsw~layer~len; 

uneigned short rasw_layer2_len; 
unsigned short lsw~layer2~lsn; 

uneigned short mavr block^len; 
unsigned short lsw^block^len; 
unsigned short nsw~msg_len; 
unsigned short lsvrjnsg~len; 
short ref; 
short dl; 



•* Data level codes to assign to output product header. 
•/ 

static unsigned short data level_16 ( NUKJ3ATAJUBVELS ) - 
< 0x8002 , OxOOOSTOxOOOa, 
OxOOOf , 0x0014 , 0x0019 , OxOOle , 0x0023 , 0x0028 , 0x002d , 
0x0032 , 0x0037 ,0x003c, 0x0041 , 0x0046 , 0x004b >; 

short pr_type ; 

unsigned short product color; 
unsigned short grid_coIor; 

short flag_off - <0>; 
short f lag~on - {1>; 

short nun; 

short prld; 

short deetlnation_id; 

short num blocks; 

short divider; 

long latitude_origin; 

long longitude origin; 

double phi center deg - {38. 0>; 

double lanoa centir_deg - {-98 ,0); 

d uble lat origin deg; 

double 1 n~origln deg; 

unsigned sH rt msw latitude; 

unsigned short law latitude; 
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unsigned short m»w_longltuda ; 

unsigned abort lsw~ longitude; 

short height_vaiue7 

short sequence^num; 

short scan_num7 

short elevetion_num; 

short eym_block~id; 

short nunflayers; 

unsigned short raater_opcode; 

unsigned short raster~constl; 

unsigned short raster~const2 > 

unsigned short raater~const3; 

char »out str home - "/horns/ - ; 
char *out~st r~site - "NATIONAL/" ; 

short num_lncluded; 
short num^excluded; 
short m ~ 

short tab included (MAX SITES); 
short tab^excludadfKAX^SITES J ; 

short year, month, day; 
long time_sec; 
long seconds; 
long cur seconds; 
short hour, minute; 
short hh,mm; 
long se ; 

struct tm JanJ)l_1970_atruct - <0>; /* Initialise all fields to 0. */ 

time t Jan_01_1970_t; 
time~t gmt~t;~ 

time t total_time> 
struct tm gmt_struct; 

short j_date; 
short cur_j_dataj 
short meg/jlate; 
long msg_tlme; 

unsigned short msw msg time; 
unsigned short lsw~msg_time; 

short scan date; 

long scan time; 

unsigned short msw sc anytime; 

unsigned short lsw~scan~timei 

short generat longdate; 

long generation time; 

unsigned short msw generat longtime; 

unsigned short Is Wegener at longtime; 

union! 

long 1 value; 
short e_vaiue(2)i 
> long2 short/ 

** Obtain a buffer into which existing products are read 

** and new products are built and then stored from- 

•/ 

but- (char •)mall c(KAZ DATA_SI«) ; 
if( buf— NULL) 
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< 

(void) printf( "error t mailoc failed for buf, errno«%d\n" , err no) ; 
return(l); 

/• 

Open the configuration file and read in the header. 

f pr» f open ( - conf ig . dat ■ , ■ rb ■ ) ; 

if ( fpr— NULL) 

i 

(void)printf ("fopen error for configuration file header, errno-%d\n" ,ermo> 
return(l); 

nread-fread(&eite file header, eizeof (site file header) , 1, fpr) ; 
if (nread— NODATA) ~ ~ 

i 

(void)printf ("fread error with aite file header, errno«%d\n" ,errno) i 
return(l); 

> 

/• 

*• Check the configuration file header for the correct keyword 
** and exit the program if the keyword doee not agree. 

if( elte_file_header. keyword I- site_key ) 

( void )printf ("Invalid configuration file keyword\n"); 
return(l) ; 

> 

/• 

*• Obtain the number of altes from the configuration file header and 
*+ use that number to determine the else of the configuration file 
** content a. Reed the configuration contents from the configuration 
** file. That provides a liet of eitee and indicates which are KCXKAO 
** and which are active (to be used in the mosaic). 
•/ 

num_sitee - aite_f ile_header.num_eiteef 

nread* f read (e it e file contents, nun sites*sizeof (site contents f ormat ) , 1, fpr) ; 
if ( nread*- NODATAj 

( 

( void ) print f ("error reading configuration file contents, errno«%d\n" , errno) ; 
return ( 1 ) $ 

> 

status - fcloee(fpr); 
if ( status J- 0 ) 

< 

( void ) print f ("error closing configuration file, errno«%d\n" , errno) ; 

> 

while(l) /* DO FORKVBR */ 

{ 

/• 

** Obtain the current time in seconds since 1/1/1970 which 
** is used to determine the age of existing products 
*• ae they are read in from the database. 

V 

Jan 0 1_19 70_st ruct . tm_year - 70; 
Jaa~01~1970~atruct.tmjBon • 0; 
Jan~0 l_1970~struct . taTmday • 0| 
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Jan_01_1970_t - mktlme (fcJan_01_1970_struct ) ; 

if (Jan_01_1970 t mm (tin* t)-l) 
{ 

( void ) print f (-error: mxtime( 1/1/70) error. \n"); 
return(l) $ 

> 

total tiroe-tlroe(NULL) ; 
grat_etruct - *gmtJjne(fitotal_time) ; 
year - gmt struct. tm year; 
if(year>99) 

< 

year— 100/ 

> 

month - gmt_struct . tm mon + 1; 
day * gmt struct . tm_m3ay / 
hour - gmE_etruct . tm_hour ; 
minute - gmt^itruct . tnajnin; 

( void )printf ( -Current <gmt) year»%d month«%d day-%d hour-%d min-%d\n w , 
year , mont h 9 day , ho u r , minute ) ; 

gmt_t « mXt ime { figmt_ struct ) ; 
if (gmt t — (time t)-l) 
{ 

(void)prlntf ( "errors mktime(gmt) error. \n*)/ 
return(l) ; 

} 

elee 

< 

time_s#c - dlfftime( gmt_t, JanJ>l_1970_t )$ 

cur_j_date - time_sec/86400; 
cur seconds - time sec%86400; 

> 



** Initialize the counters to track sites included in and excluded from the mosaic. 
•/ 

num_included -0; 
num excluded - 0; 



** Clear the grid in which the mosaic is to be built. 
*/ 

for( grld_y-KXH_Y| grld_y<ORIDJf_SXM| grid_y*+ ) 

1 for( grid x-MIN X; grid x<GRID_X_SIM; grid_x+* ) 
{ 

grid(grid_xl(grid_y] - NO_VALUEj 

> 

} 

/* 

** The following is the main loop to read in individual NIXRAD 
•* products and map the data to a large grid. Loop through 
** the sites contained in the configuration structure. 
V 

f r (j-MIW HUM; j<num sites ? 
/• 

• • Extract the site id to request from the c nfiguration file 
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** data. 
•/ 

site^id - site_file_contente( j).site_id; 



Process NBXHAD data here. Request composite reflectivity. 
prod_id * 38; 



** Find a path to the latest product of the specified type 
** for the specified site with the following. A nonzero 
•* status code indicates that no product is available 
** so jump to the next. 
•/ 

status - get_latest( site_id,prod id,prod_name); 
prod name (strlen( prod name) -1 \0 ' ; 
if (statue) 

* (void)printf (" product is unavailable. \n* ) ; 
num excluded**; 

tab~excluded(num_excluded-l] - site_id; 
goto jump here; 

> 



** Open the product file using the path to the latest. 

•/ 

fprprod-fopen < prod naae , ■ rb " ) ; 
if ( fpr_prod— NULL 7 

4 (void)printf (" Unable to open product. site-%d\n" , slte_id) ; 
nun excluded**; 

tab~excluded{num excluded-1] - aite_ld; 
pr intf ( • path«%s\n- , prodname ) ; 

goto jump here; 

> 



** Read in the product header. 
•/ 

nread- freed (buf ,header_sise, 1, fpr_prod) ; 
if< nread—NODATA ) 

1 (void)printf ("product read error, errno«%d\n%errno) ; 
nun excluded**; 

tab~excluded(nun_excluded-l) - eite_id; 
printf ( ■ path-%s\n" , prod_name ; ; 

goto jump here; 

> 

/* 

»* Bxtract data from the product header. First position the start 

• • index variable at the beginning of the buffer. 

•/ 

st_i^dex« (unsigned short *)buf; 

index-st_index ; 
mag_code» * index ; 
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Extract tha oparational moda and voluma covaraga pattarn from 
** tha product haadar. Oaa thia information to chack if tha radar 
»• vaa in a valid moda. If tha radar ia in maintananca moda 
»• than akip tha product and jump (goto) to tha and of tha loop. 
•/ 

indax-at_indax+16; 
opjabda* * indax ; 

indaV«at_lndax+17 ; 
vcp-* indax; 

if < <op_modal-l)M(opjnodai-2) ) 

* i {void)^rintt( m invalid opjnoda«%d vcp-%d\n", 
■ opjooda, vcp) / " 

num axcludad++; "~ 

tab excludadfnum axcludad-1] - aita_id; 
fer int f ( • pat h-% a \ n* , prod_nama ) ; 

foto jump hara? 

> 

/* 

*• Extract tha input product aita id numbar. Compara with tha raquaatad 
** aita id numbar. If tha axtractad and raquaatad aita id 'a do not 
** match raport an arror and jump (goto) tha naxt aita 'a product. 
•/ 

indax- at_indax*6 > 
input_aourca_id«* indax ; 

if (aita_fila_contantaljl.aita_id I- input_aourca_id) 

1 (void)printf (• raquaatad t <%d) producti<%d) aita miamatch\n - f 
aita_ld, input_aourca_id) t 
num axcludad++; ~ 

tab~axcludad[num axcludad-1) • aita^id; 
pr lntf ( " path-%a\n" , prod_nama ) j 

goto jump hara i 

> 

** Extract tha data (in daya ainca 1/1/70) and tha tima (••conda) 

" from tha product haadar. Uaa thia to datarmina tha total tima 

•* in aaconda ainca Jan 1, 1970. Thia ia uaad in a compariaon with tha 

** currant tima to datarmina tha "aga- of tha product. 

V 

indax-at lndax+20; 
j_data«*Tndax# 

indax-at indax+21; 
timal-*indax; 

indax-at lndax+22; 
tima2»* Indax; 

tot_aac-(j_data*86400)+(timal*65536)+tima2r 
if ( tot_aaot ima_aac ) 
aga minut mm - 0; 

ala 

* aga aac nda - tima a c - tot_aacf 
ag jainutaa - aga_aaconda/60j 
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> 

** Check the age of the product. If the product is too old then 
jump to the next aite. 

•/ 

if (agajninutea > max_agejninutea) 

* (void) print f (" product is too old: agejninutea»%d\n" , 

agajninutea) ; 

hum excluded**; 

tab~excludedtnum excluded-1) • aitejLd; 
printf (" path«%a\n",prod_name) ; 

goto jump_here; 

) 

/* 

** Determine the length of the data layer (in bytea). 

*/ , 

index-et lndex*66; 
lengthl-* index ; 

index-et_indax+67 ; 
length^ ■* index ; 

data_aise»( long )( (lengthl*65536) + length2 ); 
if ( deta_si»e > KAX_DAT*_SZn ) 

(void) printf ("product data site too large to read"); 
num excluded**; 

tab~excludad[num axcluded-1) - aite_id; 
printf ( ■ path-%e\n" , prod_name) ; 

goto jump here; 

> 

/• 

** Read in the reat of the product with tha following. 
•/ 

nread-f read (buf , data sixa, 1, f prjprod) ; 
if ( nread— NOOATA ) 

1 (void) printf ("error reading product file, errno-%d\n",errno) ; 
nun excluded**; 

tab~excluded(num exciuded-1) - aite_id; 
printf (" path-%e\n",prod_name) ; 

goto jump here; 

> 

atatua « fcloae(fprjrod) ; 
if ( atatua !» 0 ) 

{ (void) printf ("error cloaing input product file, ermo-%d\n",errno) ; 
return( 1) ; 

> 

*• Locate the poaltlon of the run length encoded raster data 
and extract the opcode and check it for validity. 

•/ 

at Index- ( una igned short »)buf; 
in3ex~st index; 
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opcode-* index ; 

i f ( opcode t -0xba07 ) 



S^ffiiSSS^J" 0 ' 1 inValid r " t#r ° PCOd# 10 inpUt P roduct «** \n% opcode,; 
tab_excluded{num excluded-lj - aits id; 
prtntf(- path-%e\n",prod name) J 

goto jump_here; ~ 



/• 

Extract the number of rowe from the product data layer* 

index+»9 ; 
nrowe-* index; 

/* 

Point to the first row of run length encoded raster data. 

index+«2; 
/• 

** Check for a valid number of rowe. 

if { ( nrows>MIN Jf ) fifr ( nrows<-PR_ROWs > , 
/• 

** Loop through the rowe of raster data. 

for ( row-MIH_Y ; row<nr ows ; row++ ) 

njbytes-* index ; 
index -M- ; 

cp-( char *, index; 
bin - NO_VMOT; 

/• 

** Loop through the bytes of RLS data in the current row. 
for(m-l;m<«n_bytes;m++> 

** Extract the current byte and determine the run and color. 

byte«*< °P )/ 
cp++; 

run-(( short )byte) / 16; 
col-(( short )byte) % 16; 

if (runX>) 

for(k-MlN *;k<run;k++) 

if { ( bin<PR_BINS ) && ( row<PR_ROWa > ) 

product f bin] (row) -col; 

> 

else 

{ 

/• 
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** In the unlikely •vent that the bin or row numb r 
** exceeds the array limit then report an error and 
jump to the next: site 'a product. 

*/ 

(void) print* ("error: bin«%d row-%d\n",bin,row) ; 
nua excluded**; 

tab~excluded(num excluded-1 1 ■ aits_id; 
print f [ m pat h«% • \ n " , pr od_name ) 7 

goto jump here | 

> 

bin**; 

> 

} 

> 

index- ( una igned ehort *)cpj 

> 

> 

elee 

% (void)printf ("error: product invalid number of rowa»%d« , nrows ) ; 
num excluded**; 

tab~exciuded{num excluded-1) -eite^id; 
printf (" path«%e\n w # prod w _name) ; 

goto jump here i 

> 

/* 

•* Next, open and read in the look up table file. First , create 
•* a string containing the look up table (lut) file name (which is 

is based on the eite number) and open the lut file for "read binary". 

•/ 

(void) sprint f ( lut_name, «lut/cr4_%3 . 3d. lut\0- , alte_f ile_contente [ j ) . alte_ld) ; 

fpr»fopen(lut name, •rb" ) ; 
if (fpr— NUIX)~ 

{ (void)printf (-error: error opening lut file (%•) for read, errno«%d\n w , 

lut_name , errno > ; 

num excluded**; 

tab~excludedtnum_excluded-l ] «sits_id; 
goto jump here; 

> 

** Next, read in the look up table (lut) film header from which the 

*• n umbe r of entries is extracted which is used to compute the 

*• sice of the lut file contents (lut else). Also, check the 

** keyword and the header else to verify that a valid "lut- 

** header was read in. 

*/ 

nread - f read ( *lut_header , siaeof { lut_header_f ormat ) , 1 , f pr ) / 
if( nread«0 ) 

* (void) printf ("error reading lut file, errno«%d\n", errno ) ; 
num excluded**; 

tab~excluded{num excluded-1) - site_id; 
printf (• path^eNn-fprod^name) > 

goto jump here i 

> 

num ntry • lut_header.num_entriee; 

if (Hum entry > EuTJOU) 

{ 
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num_entry - LUT_MAX; 



if( lut_headar. keyword 1- lut_key ) 

( void) print f ("invalid lut keyword - %d\n",lut header .keyword) ; 
num entry - 0; ~ 

> 

if ( lut header. header fixe l« sireof(lut header format) ) 
{ - ~ 

(void)printf ("invalid lut header size - %d\n-,lut header. header site); 
num entry « 0; ~ " 

> 

lut_size - num_entry*LUTJttDTH*sizeof (short) ; 
/• 

** Do the read to read in the lut iteelf. 

•/ 

nread - f read (lut, lut_size, 1, fpr) ; 
/• 

** Mow that the lut has been read in, cloee the file descriptor. 
•/ 

status - fclose(fpr); 
if( status I- 0 ) 

< 

(void)printf ("error* error closing look up table file, errno-%d\n",errno) ; 
return(l); 

> 

/* 

** Loop through all of the entries in the look up table 

** and nap radar data bins to grid boxes. 

*/ 

for(entry«MIN NUM; entry<num entry; entry**) 
/• 

** Extract the current radar product's bin and row coordinate. 
*/ 

bin ■ lut ( entry ) ( BIN X IDX); 
row - lut [ entry )[BnQC IDX l' 

/• 

** Check that the bin and row are valid Indlcles in the 

** product data array (that is, they are valid). 

•/ 

if((bin>»MXN X)M(bln<PR BINS ) 6ft ( row>»MIBJ r)**(row<PR ROWS)) 

/* 

** Extract the product color (data level) value from the 

** product array for the current bin and row. 

V 

product_co lor ■ product ( b in ) ( row ) j 
/• 

Extract the corresponding x and y grid coordinates. 

•/ 

grid x • lut l entry] [GRID X IDX); 
grld^y » lut [ entry ] [ GRID Jf~IDX ) ; 
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aiaa 

< 
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/• 

** Chack tha currant x and y grid coordinataa to varify that 

*• tha valuaa ara valid array indiciaa. 

•/ 

i£((grid x>»MIN X)&*<grid x<GRID X SIZB)6fi 

" <gri3 y>-MIN Yfaa ( gri3_y<GRID_Y_5 XZB ) ) 

{ 

/• 

*• Bxamina tha axiating grid box for a vaXua. 
•/ 

if< grid(grid_x}(grid_y) — NO_VALUB ) 

{ grid [grid x){grid _y) - product_color ; 
> 

aiaa 

i 

/* 

** In tha caaa with an axiating valua in tha grid tha 
*« a comparison ia naadad batvaan tha axiating grid 
•* vaXua and tha currant product coXor vaXua. 
*• Map tha product coXor whan it ia graatar than tha 
axiating grid coior vaXua. 

*/ 

grid_coior » grid(grid_x) (grid_y] j 
if <product_coior > grid_coXor) 

{ grid (grid xj[gcid_y) ■ product_eoiorj 

> 

> 

> 

aXaa 

< 

*• Thia handiaa tha unXikaXy caaa that an x,y grid coordinata 
•* waa not a vaXid array indax. For now r print tha vaXua 
** and go to tha naxt Xut tabXa antry. 
•/ 

(void)printf <-antry%d grid_x-%d grid_y*%d\n- , antry, grid_x,grid_i 



** Thia handiaa tha uniikaiy caaa that a product bin row and coXuxnn 
*• waa not a vaXid array indax. ror now, print tha vaiua 
• * and go to tha naxt Xut tabXa antry, 
•/ 

(void)printf <-antry%d bin»%d row»%d\n- , antry, bin, row) i 



> 



> 



print* (• product intagratad id-%d opa»da-%d vcp-%d\n", 
r mag_coda , op aoda , vcp ) i 

printf ( * path-%a\n* , proa^naaa) i 

nua lncludad*+i 

tab2includad[nu«_incXudad-l] - aita_iai 

jumpjiarai; /• Ju*p hara to akip a product and go on to tha naxt aita« 
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** Now that tha loop to procaaa th individual aita producta 

*» into tha grid ia complata, print out tha aitea which wara includad 

•« and tha aitaa which w«ra axcludad. 

•/ 

(void)printf ("Sitaa procaieing complatai \n m ) / 
(void)printf (" number included»%d\n" ,num_included) ; 
for(n»0; n<num_includedj n+*) 

* ( void ) pr intf ( " included«%d\n" , tab_ includad I n J ) j 

} 

< void )pr intf (■ number axcludad-%d\n B , num_excluded) ; 
for(n«0| n<num_ex eluded; n-M-) 

{ (void)printf (• excluded«%d\n" , tab_excluded[n) ) t 

> 

/* 

** Tha producta hava baan mappad to tha grid. Build a moaaic product hara. 
•/ 

prid - 5; 
/* 

Poaition tha pointar at tha beginning of tha buffar in 
** ordar to ba raady to aaalgn data to tha buffar. 
*/ 

at_index« ( una ignad ahort *)buf; 

index»e t_index+0 * 

/* 

«* Build a atrlng containing tha output filenaae. 
*/ 

( void )atrcpy( prod name,out_etr_home) i 
( void ) a treat ( prod^name, outsat r_eite ) j 

** Continua to build tha product pathnaae. Tha following builda 
** tha actual fila nama in year, month, day, hour ,minute format. 
** concatanata tha filanama to complete tha product pathnaae. 

V 

( void ) apr int f ( prod f ila , - % . 2d% . 2d% . 2d . % . 2d% . 2d\0- , 

yaar , month , day , hour , minut a ) ; 
< void ) atrcat ( prod_naae , prod_f ila > ; 

*• opan tha product pathnama for writa binary with tha 
following* 

V 

fpw-fopen(prod nama^wb"); 
if (fpw—NULL) 

* ( void ) print f<*f opan for writa arror, errno*%d\n",errno) j 
raturn(l) ; 

> 

*• Dafina tha paraaatara which apacify how tha product 

«• if t ba built* 

*/ 
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number_rowe - 896; 
number~bine « 1280; 
acale_x ■ 1; 
acale~y - 1| 

etart~row - 0; 

etart^bin • 0; 

output eource id - 10000; 

lat origin dag - phi centerjieg; 

lon~origin~deg • lam3a_center_deg; 

/* 

** Than atart aaaigning tha data to tha product header. Bvary field 
• * muat gat a value to ovarwrita any exieting data in tha buffer. 

V 

mag coda ■ prid; /* Aeeign tha maaaaga coda */ 

*in3ex * meg_code; /* which ia tha product ID. */ 

index«et index+1; 

mag date~» cur_j_date; /« Aeeign currant data in daye eince 1/1/70 

*in3ex ■ mag_date; 

mag time - cur aaconda; /* Convert second* to 2 half word*. * 

long2 abort. l_value - meg_time; 

index-at index+2; 

maw mag time - long2ehort.a value (0]; 

*in3ex - mewjaag^tlme; ~ /* Aaaign currant tima (maw). •/ 

index- at_index+3 ; 

law mag tima - long2ehort.a_value(ll i 

*ln3ex « lewjneg_tima; /• Aaaign currant tima (law). ♦/ 

maw mag lan indax«et index+4; /*Save poeition o£ tha */ 
Iaw"mag~lan;indax*at3indax+5| /• of tha maaaaga length. */ 

index-* t index+6; 

•index »~output_aource_id; /* Aaaign eourca id. */ 

index-et index* 7; 
deatinatlon id • NO_VAL0B? 

•index » deetination_id> /* Aaaign deatination id (unueed) . •/ 

index«et indexed; 
nun blocJta • 3; 

•in3ex - num_blocke; /* Aeeign tha number of blocke. •/ 

index«et_index+9 ; 

divider • -lj , _ ^ t _ . 

♦index-divider | /• Aeeign block divider. •/ 

latitude origin - ( long) (lat_orlgindag* 1000. ) ; 
long2ehort.l_value - latltude_origin; 

index-it index+10; 

maw_latitude « long2ehort.e_value[0J j 

♦index « mew_latitudai /* Aaaign origin latitude (maw). */ 
index-et index-fill 

lew latitude - 1 ong2 abort. e_value(l) ; 

•ln3ex - lew_latitudai /♦ Aaaign origin latitude (law). */ 

1 ngituda rigin • (long) ( lon_origin_deg* 1000. ) i 
1 ng2ehort.l_valu • 1 ngitude_origini 

index-at index+12; 
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maw longitude - long2 abort .a_valua(0] ; 

•in3ax • maw_longltuda; /• Aaaign origin longituda (maw), */ 
indax»at indax+13; 

law longituda * long2ahort.s valua(l); 

*in3ax » law_longituda; /* Aaaign origin longituda (law). */ 

indax-at_indax+14 ; 
haight valua - NO VALUE; 

*indax~- haight_vi*lua; /* Aaaign tha haight valua (unuaad).*/ 
indax*at indax+15; 

♦indax »~prld; /• Aaaign tha product ID. */ 

indax»at_indax+16 ; 
op rooda « 2; 

•indax • op_moda; /* Aaaign tha oparatlonal moda valua. •/ 

indax»at_indax+17 ; 
vcp • NO~VALUB ; 

♦indax «~vcp; /* Aaaign tha voluma covaraga pattarn. */ 
indax»at_indax+18 ; 

aaquanca~num ■ NO VALUB; . 
♦indax »~aaquanca~num; /♦ Aaaign tha aaquanca numbar (unuaad). */ 

indax*at_indax*19 ; 

acan num"V NO VALUE; Jk m . 

♦indax « acan~num; /* Aaaign tha voluma scan numbar (unuaad). ♦/ 

indax»at indax+20; 
acan data - cur j data; 

* indax - acan_data; /* Aaaign tha acan data. •/ 

acan tima - cur aaconds; 

long5ahort.l_vaIua - acan_tima; 

indax»at indax+21; 

maw acan~tima - long2ahort.a_valua(0) ; 

♦in3ax -~maw_acan_tima; /* Aaaign tha acan tima (maw). ♦/ 
lndax«at lndax+22; 

law acan~tima - long2ahort.a_valua{l) ; 

♦in3ax •"law^acaB^tima; /* Aaaign tha acan tima (law). */ 

indax«at indax+23; 
ganaratlon data - cur_j_data; 

* indax - ganaration_data; /* Aaaign tha ganaratlon data. */ 

ganaratlon tima - cur_aaconda; 
long2ahort.i_valua « ganaration^tima; 

indax»at indax+24; 

maw ganaratlon tima * long2ahort.a_valua[0J ; 
*ln3ax - maw ganaration_timm# /* Aaaign tha ganaratlon */ 
— /* tima (maw).*/ 



indax-at indax+25; 

law ganaratlon tima - long2 abort. a^vaiuaf 1 J ; 
♦lnHax - law_ganaratlonJ;ima; /• Aaaign tha ganaratlon */ 

^x'^S^raJ /• Product Dapandant (unuaad). •/ 

iKV^S5£&! /# product 
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index»st_index+28 ; 
elevation num - NO VALOT; 

♦Index * elevstion~nuin;/* Assign the elevation number (unused). •/ 



index-et index+29; 
•index «~NO_VALUB> 

index-st index+30* 



/* Product Dependent (unused). •/ 
/• Save position of data levels. •/ 



•* Assign the appropriate set of data levels values. 



for( dl-HIN_VAtf7E; di<NUM_DATA_LEVlLS ; di++ ) 

•index - data level 16 [dl]; 
index** j 



index-st index+46r 
♦index -""NO^VALUB; 

index-st lndex+47; 
•index -~NO_VALOTi 

index-st index+48; 
•index «~NO_VALU*/ 

index-st index+49; 
•index -~NO_VAI,UK/ 

index-st index+SOj 
•index «~NO_VAL0*j 

index-st index+51; 
•index -~NO_VALUB| 

index-st_index+52 ; 
•index -~number_blns % 

index-st index+53j 
•index -~NO_VALUK > 

index-et index+54; 
•index -~NO_VALOT; 

index-st_index+55 $ 
•index —""60 # 

index-st index*56j 
•index -~NO_VALUKj 

index-st lndex+57; 
•index -~NO_VAI*OT* 

index-st index+Sfl; 
•index «~NO_VALOT> 

index-et index*59r 
•index -~N0_VALOTj 

index-st index+60; 
divld r - -1; 
•ind x-divlderi 

index • st lndex+61; 



/* Product Dependent (unused). 

/* Product Dependent (unused). 

/« Product Dependent (unused). 

/* Product Dependent (unused)* 

/• Product Dependent (unused). 

/• Product Dependent (unused) . 

/• Number of columns (bins). •/ 

/• Onused. */ 

/* Offset to symbology (msw) . •/ 

/* Offset to symbology (lew). */ 

/♦ Offset to graphic (msw). •/ 

/• Offset to graphic (law). •/ 

/• Offset to tabular (msw), •/ 

/• Offset to tabular (lsw). •/ 

/• Beginning of symbology block •/ 

/• Assign block divider. •/ 
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•ym block id • 1; , _ . _ . 

•inoax-aym_block_id; /* Aaaign ayrabology block 10. •/ 

indax- at indax+62; 

maw_block~len_indax - index; /• Sava location of block aiae valua*/ 

indax- at indax+63; 
lew_biock~lan_indax - indax; 

indax- at_index+64; 
nun layar a - 1; 

*in3ax - numJLayara; /* Aaaign numbar o£ layara */ 

indax- at index+65? ^ m 

•indax - Hlvidar; /• Aaaign layar dividar. •/ 

indax- at indax+66; m _ * 

maw_layar~lan_lndax - indax; /* Sava location of layar langth. */ 

indax- *t_indax+67; 
iaw_layar_lan_indax - indax; 

indax- at indax+68; 
raatar opcoda - 0xba07; 

• indax"- raatar_opcoda; /• Aaaign Raatar data format opcoda. •/ 

indax- at indax+69; 
raatar conatl - 0x8000; 

*indax~- raatar_conatl; /* Aaaign Raatar data format conatant. •/ 
indax* at indax+70; 

raatar conat2 • OxOOcO; _ 

• indax"- raatar_conat2; /• Aaaign Raatar data format conatant. V 

•indax ^NO^VALOl; ' /* Aaaign placaholdar for raatar x origin. •/ 

• indax I t ^^o"vA^OT? , /• Aaaign placaholdar for raatar y origin. •/ 

indax- at indax+73; ^ i fc m/ 

•indax - T; /* Aaaign acala x (intagar portion). •/ 

•lnd« VvtoWLmV /• Aaaign acala x (fractional portion). •/ 

•indax I*?)***** 7 /*' Aaaign acala y (intagar portion). •/ 

SXS VltoySnV /• Aaaign acala y (fractional portion). V 
indax- at indax+77; 

• indax - numbar_rowa; /• Aaaign numbar of rowa. •/ 
indax- at_indax+78; 

• ln^--°rtatar jonat3; /• Aaaign packing daacriptor. •/ 

n^rindax-indax, /•Sava pojition of tha -numbar of •/ 
* - /• bytaa thia row" for tha lat row, »/ 

laymr_langth « 22; 

• • Tha following run langth ancodaa tha data in tha grid. 
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Prmv col— l; 
run-MIN VALUE; 
nbytea.MIN VALUE ; 
ind«x++/ 

cp-(char *) index; 



> 



;; loop t ^,, .u o t ts . bl „. ta th . _ 

^•""-"".^olo.oo^j.t,.,^, 

;; Mt. .„ th . t only jrr^ri^srt£i slst 

col - grid (bin J (row); 

{ _ l)«(colj. pr .^ col) , i i (run .. 15) j 

S^olcr^h.^. " 5 h# ginning o£ . row wh .„ 

(which u srsLnS d jr h * n ~ r ~ n * run\: h « 

- combination i.^-T^JT 1 ? • Th. run/color 

position 



cp++; 

run-i; 

pr«v_col - col; 

•In 

i 



/* Otherwise, update i„ .-i 

« 9rid box UifZltZ.™ • 3,i « ti «9 «n and go on to th. n .xt 



run**-; 

Pr«v_col - col; 



i* (run>l> 

{ 



- St r. th. l„t -run color . coabln . tion ^ ^ 
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** utput buffar with th following* 
•/ 

byta» ( run* 16 ) +prav_col j 

*(cp)»byta; 

cp++; 

nbytaa-M-/ 

> 



** If tha number of bytaa ia odd than pad 
** ona mora byta with a valua of taro. 

if(nbytaa % 2) 

{ 

byta-MIN VALUE; 
*(cp)«byta; 
cp++; 
nbytaa++; 

> 



** Aaaign tha numbar of bytaa for thia row to tha 

** buffar in tha propar location. 

•/ 

* (nbytaa indax > -nbytaa j 
layar_langth+»nbytaa ; 



** Updata layar langth to includa 2 bytaa 
** for tha nbytaa valua. 

*/ 

layar_langth+-2 1 

/• 

•* Advanca tha indax point ar ao that it point a to tha 
location of tha naxt row* a byta langth. 

•/ 

indax- ( una ignad ahort *)cp# 
nbytaa indax- indax; 

> 

** Now that all of tha rowa in tha grid hava baan procaaaad, 

** convart tha layar langth valua to 2 -halfworda- 

** and aaaign thoaa to tha layar alaa poaition in tha 

•« product aymbology layar ona haadar. 

•/ 

long2ahort.i valua • layar_langth; 
maw layar lan - long2ahort.a valua(O); 
l*v~layar~lan » long2ahort.a_valua(l)i 
•maw layar lan indax - maw_layar_lan; 
* 1 aw~l ayar^lan^indax • law_layar_lani 

block_langth - layar_langth+16 ; 

** Convart tha block langth valua to 2 "halfworda* 

*• and aaaign thoaa t tha block langth poalti n in tha 

•* product aymbology block haadar. 

•/ 
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long2ah rt.l value - block_l ngth; 
maw block len « long2 short . e_value J 0 J ; 
le%Tblock~len - iong2ahort.s value(l]; 
•maw block" Ion indax - mew_biock_leni 
*lsw~block~len~lndex - iaw_blockJLem 

/• 

** Convert tha meeaage length value to 2 -halfworda- 
*• and aaalgn thoae to the meaaage length position xn the 

* * product message header block. 
•/ 

total length - block_length+18+102; 
meeeage_length ■ total_lengthj 
long2ehort.l valus - meeaage_lengthj 
maw msg len » long2ahort. a_value(0] ; 
ls%/"msg~len - long2ahort.a_valuetl]; 
•mew mag len index - mewjnag_len; 

* law^mag^len^index • lawjnag_lsm 

/* 

** Writs ths product to ths fils and cioas ths fils. 
*/ 

nwrits - fwrlts(buf ,total_length,l, fpw) j 
if (nwrito—NODATA) 

* ( void )printf ("error writing product file, errno«%d\n- ,errno) j 
rsturn [1) i 

statue » fflush(fpw)j 
It { status l« 0 ) 

* (void )printf ("error: error flushing product fils, errno»%d\n- ,errno) $ 
return (1) j 

statue - fcloae(fpw); 
if ( status !• 0 > 

{ (void)printf ( -error t error cioaing product file, errno«%d\n- ,errno) t 
return < 1) / 

> 

print f (-Product is complsts. \n\n-); 
( void ) f f lush < st dout ) ; 
(void) f flush ( stdsrr ) ; 

rsturn (0) ; 

* Routins to obtain latsst vsrsion of ths requestsd product for ths requested aite. 
[it get_lateet< short sits^id, ahort prod_id, char «s> 

* This passes a string containing the cc»mplete pathnaas 

* to thi lateet product in the data baae °*^******" tfimA 

* by the product id and for the sits spsclfisd by sits id. 

* A statu, cods is rsturnsd by the routine, ^«?^°g r , thm 

* return valuesi rsturn • 0 -> product name fDun f n ( ^LS^ 
* return - 1 -> no product found in directory. 

* return * 2 -> invalid product id number. 

* return • 3 -> invalid aite id number. 

* return - 4 -> popen failed. 
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V 

char •teat; 

char commandt512l; 

char »t ring [512], 

conat char *commandl « « /u«wr< « /i 

FIL* *cfd; 
ahort j; 

» copy th. f ir . t part of th . cofflmand tQ ^ 
( void ) .trcpy ( coonand , commandl ) ; 

SSETSiJ? (h °"-' di "<*°*V .trin, to th. 

(votd).trcmt(conmand,.tring hoa.), 
/• 

s.« ch for th . .p. clfied , lt . td with th# foiiewia9 ieop 

for(j-0/ J<KDK_SITK_PATH/ 

it(.it._p at h[j)..tt._id — .it._id , 

- jump ?br!j? r S 3 £ rooT™ 1 " trin9 Md th « n 

> 

/* 

•trinofO)-'\0'; 

Jl° id) ** re Py ( •tring) | 
r.turn(2); 

ju«pl_h.r.,, /. juap to h . r . Md contlnu#< #/ 
/• 

ISoffid'Td 19 •• areh " th « of product id', for th 

J r O-°> J<H0M_PROD_PATH; 
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if (prodjpathtj ).prod_id -« prod_id ) 
/• 

** When the product id la found, concatenate tha product 
** nana path »tring to tha command atring and than 
** jump (break) out of tha look. 



( void ) at rcat < command , prod_path ( j ] . prod text ) ; 
goto jump2_herej 

> 

/• 

** The following ia executed only when the product id ia 
** not found. Null the output atring, copy it to the output, 
** and return from the routine with a nonaero value indicating 
** an invalid path. 



stringfO]»'\0' ; 

(void) etrcpy< a, atring) ; 

return(3); 

jump2_hereif /* jump to here and continue. »/ 
/* 

** Concatenate the laat part of the command on to the end of the 
** command atring. 

V 

( void) at rcat( command , command*) ; 
/• 

** open a pipe and execute the command under the "eh- ahell. 

** The reault la returned in a atring pointed to by the pointer "cfd*. 



cfd~popenf command, "r") ; 
if (Cfd— KUIX) 

< 

/* 

** Thia handlea the unlikely caee that the popen function failed. 

** Print a meaaage and return with a nonaero atatue code* 

•/ 

( void) print f ("get lateat error -> cfd ia (TOLL: errno-%d\n" ,«rrno) ; 
return (4); 

> 

/* 

** Get the atring with the following. Remember that fgete returna 
** a atring which include a newline character at the end (which muat 
•* be removed by replacing it with a KUL character. 
•/ 

etring[0]«'\0' $ 

( void ) a t r cpy ( a , atring ) ; 

teat- f get a ( atring, eiaeof ( atring) , cfd) ; 

pcloae(cfd) ; 

/• 

** Check the fir at character of the atring for a null. 
•/ 

if (etringlO) — '\0' ) 
{ 
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/• 

* # This handle* a string with the first character a null. 
"It Man* that the requested product was not found. 
+* Tha null string ia paaaad and tha routine returns with a 
** nonzero value. 

V 

( void ) strcpy ( a , string ) ; 
return* 1) ; 

> 

/• 

** When the string is valid, replace the laat character (a carriage return) 

with the null character. This insures that just the path name ia 
** returned to the calling routine. 

string ( strlen ( string ) -1 } - ' \0 ' i 
/* 

** Copy the string to the output string which is paased to the calling 
*• routine. The value of zero is returned with return indicating 
** that a valid product waa found. 
•/ 

( void ) atrcpy < a , atr ing ) ; 

return(O) ; 

> 

The apparatus and method of the present invention 
converts a radar data bin (or the like) to a grid location 
or box in a resulting mosaic, rather than vice versa* If 
the reverse is done, occasional radar data bins could be 
missed. 

The present invention can construct a mosaic using any 
data that is geographical or otherwise territorial (in its 
broadest sense) in nature or relating to merging into one 
mosaic points from a plurality of surfaces. However, it is 
strongly preferred that such data be in a raster format. 
Data from different sites can thereby be combined to 
produce a single mosaic. 

The present invention can be practiced using prebuilt 
lookup tables, or by performing the reprojection 
calculations in real time. However, use of lookup tables 
is preferred, for faster operation and quicker response in 
real time. 

Some of the many advantages of the invention should 
now be readily apparent. For example, apparatus and method 
are provided for combining of data from multipl sources 
int one r m re m saics such as would be sultabl for 
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utilization by nd users. Such apparatus and meth d is 
capable of accurately and efficiently constructing such 
mosaics. Also, the location of any holes in such a mosaic 
can thereby be identified. The present invention also 
provides for filling such holes. 

Thus there has been provided apparatus and method for 
combining data from a plurality of sites into a single 
image covering a desired geographical area* Such apparatus 
and method are capable of combining, in real time, radar 
data (such as weather radar data) or other data in digital 
format from multiple radar sites into a mosaic covering a 
geographical area. Such apparatus and method are 
furthermore capable of combining data from a plurality of 
sources into a single image that can be run on a general 
purpose computer such as a personal computer. Such 
apparatus and method are capable of producing a mosaic of 
weather radar data that can provide location and intensity 
of precipitation and other meteorological phenomena 
throughout a large geographical area. When lookup tables 
are used, the present invention can accurately combine a 
substantial amount of data from a plurality of desired 
sites fast enough to provide a current depiction of the 
weather or other phenomena. 

Obviously, many modifications and variations of the 
present invention are possible in light of the above 
teachings. It is therefore to be understood that the 
foregoing embodiments are presented by way of example only 
and that, within the scope of the appended claims and 
equivalents thereto, the invention may be practiced 
otherwise than as specifically described. 
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What is claimed is: 

CIAIMB 

1. A method for inserting a first set of data having a 
first coordinate system into a data mosaic having a second 
coordinate system which may be different from the first 
coordinate system, comprising the steps of: 

converting the location of one item of the first set 
of data into a first distance from a first reference point 
in the first coordinate system; 

converting the first distance to a location in a third 
coordinate system, wherein the location of the first 
reference point in the third coordinate system is known; 

converting the location in the third coordinate system 
to a second distance from a second reference point in the 
second coordinate system, wherein the location of the 
second reference point in the third coordinate system is 
known ; and 

converting the second distance to a location in the 
second coordinate system. 

2. A method for inserting a first set of data having a 
first coordinate system into a data mosaic having a second 
coordinate system which may be different from the first 
coordinate system, comprising the steps of: 

forming for the first set of data a table having a 
plurality of entries, each such entry having a location in 
the first coordinate system having any of the first set of 
data, and each location in the second coordinate system 
corresponding to that location in the first coordinate 
system, said forming step comprising the steps of: 

converting the location of one item of the first set 
of data into a first distance from a first reference point 
in the first coordinate system; 

converting the first distanc t a location in a third 
coordinate system, wh r in the locati n of the first 
r ferenc point in the third coordinat system is known; 
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converting the location in the third coordinate system 
to a second distance from a second reference point in the 
second coordinate system, wherein the location of the 
second reference point in the third coordinate system is 
known; 

converting the second distance to a location in the 
second coordinate system; and 

after said forming step, using the table, converting 
each of a plurality of locations in the first coordinate 
system having at least some of the first set of data to one 
or more locations in the second coordinate system. 

3. A method for inserting a first set of data having a 
first coordinate system into a data mosaic having a second 
coordinate system which may be different from the first 
coordinate system, comprising the steps of: 

converting first locations for the first set of data 
in the first coordinate system into second corresponding 
locations in the second coordinate system; 

defining an array of data of a first value in th 
second coordinate system, the array including all of the 
second locations; 

changing the value of each of a first subarray of the 
array of data from the first value to a second value, the 
subarray consisting of all locations in the second 
coordinate system that could be mapped into any of the 
first locations in the first coordinate system; 

changing the value of each of the second locations in 
the array of data from the second value to a third value; 

for a location in the array having the second value, 
determining the positionally closest corresponding one of 
the first locations; and 

inserting the value of the positionally closest 
corresponding one of the first locations in a location of 
the second coordinate system corresponding to the location 
in the array having th second value. 



WO 96/15504 



PCT/US95/14479 



- 73 - 

4. Apparatus for ins rting a first s t of data having a 
first coordinate system into a data mosaic having a second 
coordinate system which may be different from the first 
coordinate system, comprising: 

means for converting the location of one item of the 
first set of data into a first distance from a first 
reference point in the first coordinate system; 

means for converting the first distance to a location 
in a third coordinate system, wherein the location of the 
first reference point in the third coordinate system is 
known ; 

means for converting the location in the third 
coordinate system to a second distance from a second 
reference point in the second coordinate system, wherein 
the location of the second reference point in the third 
coordinate system is known; and 

means for converting the second distance to a location 
in the second coordinate system. 

5. Apparatus for inserting a first set of data having a 
first coordinate system into a data mosaic having a second 
coordinate system which may be different from the first 
coordinate system, comprising: 

means for forming for the first set of data a table 
having a plurality of entries, each such entry having a 
location in the first coordinate system having any of the 
first set of data, and each location in the second 
coordinate system corresponding to that location in the 
first coordinate system, said means for forming comprising: 

means for converting the location of one item of 
the first set of data into a first distance from a first 
reference point in the first coordinate system; 

means for converting the first distance to a 
location in a third coordinate system, wherein the location 
of the first reference point in the third coordinate system 
is known; 
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means for converting the location in the third 
system to a second distance from a second reference point 
in the second coordinate system, wherein the location of 
the second reference point in the third coordinate system 
is known; 

means for converting the second distance to a 
location in the second coordinate system; and 

means, using the table, for converting each of a 
plurality of locations in the first coordinate system 
having at least some of the first set of data to one or 
more locations in the second coordinate system. 

6. Apparatus for inserting a first set of data having a 
first coordinate system into a data mosaic having a second 
coordinate system which may be different from the first 
coordinate system, comprising: 

means for converting first locations for the first set 
of data in the first coordinate system into second 
corresponding locations in the second coordinate system; 

means for defining an array of data of a first value 
in the second coordinate system, the array including all of 
the second locations; 

means for changing the value of each of a first 
subarray of the array of data from the first value to a 
second value, the subarray consisting of all locations in 
the second coordinate system that could be mapped into any 
of the first locations in the first coordinate system; 

means for changing the value of each of the second 
locations in the array of data from the second value to a 
third value; 

means, for a location in the array having the second 
value, for determining the positionally closest 
corresponding one of the first locations; and 

means for inserting the value of the positionally 
closest corresponding one of the first locations in a 
location of the s cond coordinate syst m corresponding to 
the location in the array having the second value. 
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